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SECTION  1 


INTRODUCTION  AND  SUMMARY 


SUMMARY 

This  report  presents  results  from  a  Naval  Research  Laboratory-sponsored 
investigation  of  dual  polarization  gyros.  The  dual  polarization  concept  is 
a  new  fiber  optic  gyro  configuration  designed  to  provide  independence  from 
principal  gyro  bias  error  sources.  Experimental  results  indicate  an  order 
of  magnitude  improvement  in  bias  error  performance. 

Important  results  from  a  new  modeling  capability  are  also  included.  In 
addition  to  examination  of  dual  polarization  gyro  performance,  a  number  of 
fiber  gyro  configurations  have  been  examined.  Results  of  this  effort 
include  development  of  an  improved  understanding  of  fiber  optic  depolarizers 
and  identification  of  a  previously  unreported  error  term  for  gyros  using 
3  dB  distributed  couplers. 

GYRO  ERROR  TERMS 

The  most  significant  difference  between  rotation  sensors  and  other  sensors 
such  as  hydrophones  or  magnetometers  is  the  requirement  for  absolute 
measurements.  The  hydrophone,  for  example,  only  needs  stability  for  periods 
of  time  on  the  order  of  a  fraction  of  a  second.  Gyros  must  maintain  an 
absolute  reference  from  hours  to  years,  depending  on  the  mechanization. 

This  absolute  reference  requirement  is  easier  to  meet  with  the  Sagnac 
interferometer  currently  used  for  the  active  ring  laser  gyro  than,  for 
example,  a  Mach-Zehnder  interferometer.  Both  geometric  paths  in  the  Sagnac 
are  identical.  Components  in  the  Sagnac  interferometer  are  sensitive  to 


environmental  effects,  giving  rise  to  non-reciprocal  phase  shifts  and, 
hence,  a  non-ideal  path  for  light  traveling  in  the  two  directions.  This 
reduces  the  stability  of  the  device,  causing  gyro  bias  errors. 

It  is  important  to  note  that  despite  the  popularity  of  shot-noise-based 
error  calculations,  the  principal  problems  in  gyro  performance  have  always 
been  component-related.  One  of  the  most  significant  components  in  terms  of 
stability  is  the  phase  modulator,  which  is  used  to  achieve  a  V2  offset 
between  the  two  counter  rotating  beams.  The  environmental  sensitivity  of 
phase  modulator  components  is  especially  important  in  the  military  gyro 
environment,  where  one  normally  requires  instant-on  performance  in  a  rough 
environment  without  temperature  compensation.  Honeywell's  effort  has  been 
directed  at  this  critical  phase  modulator  element.  Evaluating  a  new  dual 
polarization  concept  for  reducing  phase  modulator  errors  is  the  basis  of  the 
program. 

EXPERIMENTAL  MEASUREMENTS 

Tests  reported  in  this  docunent  were  performed  with  a  rate  table  which 
supported  the  optical  experiment  on  a  small  honeycomb-reinforced 
ferromagnetic  surface.  All  the  power  and  signal  connections  to  the  table 
were  passed  througn  slip  rings,  allowing  for  complete  freedom  for  rotational 
measurements  over  a  ^120  deg/sec  range.  The  flexibility  of  this  test  set-up 
allowed  convenient  interchanging  of  parts  and  reconfiguration  of  the  gyro. 

A  dedicated  computer  facility  was  used  for  experimental  measurements.  The 
computer  provided  real  time  storage  of  data  on  disk  files.  This  computer 
was  configured  to  provide  a  display  of  gyro  statistics  during  drift  runs, 
accounting  for  necessary  data  normalization  and  processing.  This  greatly 
assisted  the  experimental  set-up  by  providing  quantitative  data  on  gyro 
performance.  Both  analog  and  computer  driven  x-y  plotters  were  installed  as 
part  of  the  measurement  instrumentation  package.  Plots  of  gyro  performance 
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contained  in  Section  4  of  this  report  were  generated  directly  by  the 
dedicated  computer  facility  used  for  the  experiments.  Software  listings  for 
the  data  processing  and  display  functions  are  included  in  Appendix  A. 

In  addition  to  the  experimental  efficiency  provided  by  this  computer  set-up, 
a  number  of  unique  diagnostic  capabilities  were  available.  In  one  case,  a 
laser  noise  term  was  identified  by  Fourier  analysis  of  laser  output  and  gyro 
detector  signals.  Because  a  large  number  of  parameters  could  be  recorded, 
it  was  possible  to  correlate  variables  such  as  temperature  with  other 
experimental  data.  This  capability  was  used  to  provide  experimental 
diagnostics  as  well  as  output  plots. 

The  heart  of  the  issue  investigated  in  this  contract  effort  was  the 
feasibility  of  a  dual  polarization  gyro  concept.  The  program  objective  was 
to  demonstrate  that  the  dual  polarization  gyro  could  provide  some  degree  of 
independence  from  bias  error  problems.  The  greatest  difficulty  in 
maintaining  bias  stability  was  expected  to  be  the  temperature  dependence  of 
the  phase  modulator.  The  experimental  program  described  in  this  report  had 
a  goal  of  evaluating  the  dual  polarization  gyro  performance  in  terms  of 
these  temperature  problems. 

Gyro  operation  was  examined  as  a  function  of  temperature.  The  single 
(conventional)  polarization  performance  and  dual  polarization  performance 
were  compared.  Data  presented  in  Section  4  of  this  report  indicate  that 
with  temperature  variation,  the  dual  polarization  configuration  has  an  order 
of  magnitude  improved  bias  stability. 
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FIBER  SENSOR  MODELING 


A  numerical  solution  for  gyro  performance  was  completed  using  a  Honeywell 
developed  interferometric  model.  The  unique  properties  of  this  simulation 
are  its  applicability  to  arbitrarily  large  and  complex  systems,  including 
dual  polarization,  multiple  paths,  and  many  optical  components.  The  optical 
system  to  be  modeled  is  described  by  input  lines  at  execution  time,  rather 
than  being  part  of  program  logic.  This  "compiler''  approach  allows 
flexibility  in  the  modeling  of  complex  optical  systems. 

A  Jones  vector -matrix  representation  was  used  for  the  optical  components 
being  analyzed.  Wave  amplitudes  were  represented  by  complex  two-vectors  for 
both  polarizations.  As  the  interferometric  system  is  modeled,  the  total 
wave  amplitude  and  absolute  phase  along  each  path  are  retained  in  a  complex 
propagation  factor.  The  Jones  vector  representation  is  not  normalized  in 
the  conventional  manner,  allowing  preservation  of  absolute  amplitude  along 
each  path  of  the  model. 

In  addition  to  the  dual  polarization  gyro,  a  number  of  important  optical 
components  have  been  analyzed;  results  are  reported  in  Section  3.  They 
include  a  narrow-band  birefringent  depolarizer,  birefringent  fiber  sections, 
and  a  simple  distributed  coupler  gyro. 

The  work  with  the  depolarizer  has  led  to  important  understanding  of  how  this 
component  can  be  designed  for  gyro  systems.  Important  design  parameters 
such  as  the  characteristics  of  depolarizer  section  length  have  been  examined. 

The  effect  of  imperfect  splicing  of  fibers  has  been  considered  both  in  the 
case  of  the  dual  polarization  gyro  and  as  a  generic  gyro  component.  The 
results  of  this  study  indicate  that  birefringent  fibers  must  be  carefully 
aligned  to  preserve  desired  gyro  properties. 
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A  fourth  example  which  was  modeled  is  the  simple  gyro  configuration,  which 
includes  only  a  distributed  coupler  and  sensing  fiber.  In  this  case,  a 
significant  result  was  the  discovery  of  the  large  error  term  coming  from 
distributed  coupler  temperature  sensitivity.  For  small  fiber  optic  gyros, 
this  error  can  result  in  drift  terms  which  amount  to  hundreds  or  thousands 
of  degrees  per  hour  for  temperature  changes  on  the  order  of  10  degrees 
Celsius. 
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SECTION  2 


THEORETICAL  DISCUSSION 


In  order  to  achieve  fundamental  performance  limits  in  a  fiber  optic  rate 
sensor,  thermal  mechanical  and  magnetic-induced  errors  must  be  removed  from 
the  device.  Components,  including  the  source,  detectors,  couplers  and 
modulators  (both  phase  and  frequency),  represent  one  class  of  error  terms 
that  play  a  major  role  in  reducing  the  performance  of  fiber  optic  rate 
sensors.  These  error  sources  are  dependent  on  the  measurement  concept  and 
configuration. 

A  second  more  fundamental  class  of  errors  exists  in  the  fiber  itself.  The 
non-reciprocal,  scattering,  and  birefringent  properties  of  the  fiber  are 
truly  the  practical  limiting  error  sources  in  the  fiber  optic  rate  sensors. 
In  order  to  achieve  the  projected  noise  performance  limits,  the  mechanical, 
thermal,  and  magnetic  errors  of  the  fiber  must  be  overcome. 

A  third  problem  must  also  be  resolved  to  achieve  the  best  noise-limited 
performance:  the  electronic  readout  technique  used  to  detect  fringe  motion 
must  be  capable  of  detecting  shifts  in  the  range  of  10  ^  to  10  ® 
radians.  This  represents  a  significant  technical  challenge  that  must  be 
faced  in  all  interferometric  devices. 

In  the  following  subsections,  a  generalized  development  of  the  Sagnac 
interferometer  output  is  presented,  allowing  consideration  of 
component-induced  bias  errors.  An  alternative  measurement  scheme  involving 
two  polarizations  is  discussed  in  detail.  This  dual  polarization 
configuration  has  been  demonstrated  experimentally.  Results  are  presented 
in  Section  5  of  this  report. 
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FIBER  OPTIC  RATE  SENSOR  THEORY 


Fiber  optic  rate  sensing  is  based  on  the  Sagnac  interferometer.  The  classic 
treatment  of  the  Sagnac  effect  by  E.J.  Post1  shows  that  the  phase 
separation  between  the  two  counter-rotating  beams  is  proportional  to  the 
input  rotation  rate  and  is  of  the  form 

A<t>  =  a 

XC 

where 

K  =  number  of  turns  of  fiber 

A  =  area  enclosed  by  a  single  fiber  loop 

X  =  free  space  wavelength  of  the  light  source 

c  =  velocity  of  licfit  in  free  space 

fl  =  input  rotation  rate,  normal  to  the  loop  plane. 

This  is  the  classical  Sagnac  expression.  In  order  to  better  understand  the 
effects  of  errors  induced  in  the  device,  it  is  necessary  to  develop  a  more 
general  form  of  the  Sagnac  expression  that  will  allow  for  explicit  examination 
of  the  error  terms. 

The  primary  component  of  a  fiber  optic  rate  sensor  is  the  closed  light  path 
that  is  formed  into  an  interferometer.  Light  moving  along  the  path  will 
experience  a  phase  change  proportional  to  the  angular  rate  of  the  path  in 
inertial  space.  The  sense  of  the  phase  change  is  determined  by  the  direction 
of  the  light  motion  with  respect  to  the  angular  rate.  Therefore,  if  a  light 
beam  is  split  into  two  beams  that  travel  in  opposite  directions  along  the 
light  path,  these  two  beams  will  experience  equal  and  opposite  phase  changs. 
The  phase  changes  are  also  proportional  to  the  projected  enclosed  area  whose 
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normal  is  parallel  to  the  axis  of  rotation.  As  a  result,  the  axis  of  rotation 
can  be  determined  by  using  three  independent,  mutually  orthogonal,  closed 
light  paths. 

Consider  two  interferometer  beams  that  have  plane  wave  phase  fronts  in  free 
space  and  in  the  interferometer  path  medium.  The  electric  fields  are 
described  as: 

E1  =  E01  Exp  i^kl  *  r  “  “i  t  + 

E2  =  EQ2  Exp  i(k2  •  r  -  oj2  +  m2) 

where 

=  wave  vector  at  frequency 
=  arbitrary  phase 

Eq^  =  electric  field  amplitude/polarization  vector 
"f  =  displacement  along  fiber. 

The  accunulated  phase  expression  for  the  light  beam  after  circulation  around 
the  path  is: 

x 

<t>£  =  •  "dr  -  /a)  dt  +  m1  (i  »  1,  2) 

o  i 

where  t  is  the  time  required  to  traverse  the  path.  For  a  path  length 
L,t  =  nL/c,  where  n  =  the  index  of  the  medium  and  c  =  the  velocity  of 
light  in  space.  For  a  system  in  which  the  observer  is  co-moving  with  the 
interferometer,  the  spatial  integral  is  a  constant  and  only  the  temporal 
integral  is  of  interest.  When  the  interferometer  path  is  subjected  to  a 
rotation  generated  by  velocity  vector  V,  the  change  in  the  light  beam 
transit  time  can  be  written  as* 
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tj  =  —n  /nt  (1  -  a)  V  •  dr  (i  =  1,  2) 
1  (T  1 


where  a  is  a  drag  coefficient  similar  to  the  Fresnel-Fizeau  coefficient. 

The  relative  phase  in  time  between  the  two  counter-directional  beams  can  be 
written  as 

Ti  t2 

A<J>(t )  =  ^(t)  -  <J>2(t )  =  f  oijdt  -  S  oijdt  +  (m^  -  m2) 


where  x^  =  x^  +  x^  +  t.  If  we  assume  that  and  are 
constant,  then  with  the  refractive  index  and  drag  coefficient  uniform 
throughout  each  path, 


L  LV  I  2 

=  —  4>2n2 )  +  2  **h(l  —  o^) 

+  n|  u>2  (1  -  o^)  j  +  (m1  -  (i>2 ) t  +  (m^  +  m2). 


Consider  now  that  the  path  is  a  circle  of  radius  R  with  N  turns  of  fiber 
around  it.  The  tangential  velocity  V  can  be  expressed  in  terms  of  an  angular 
rate  as  v  -  R n,  and  the  path  length  is  L  =  2ttRN.  Using  these  expressions 

O 

and  assuming  the  drag  coefficient  to  be  of  the  form  a  =  (1  -  n  ),  the 
phase  can  be  written  as 

=  2t^2r2  (c^  +  o^)  +  (u^  ^  rij) 

+  (cu^  -  w2)t  +  (m^  -  m2). 

points  P  and  Q.  Different  choices  of  parameters  for  best  distribution  on  the 
sphere  are  very  easily  explored  using  the  computer  program.  Design  iterations 
are  efficient  because  of  the  flexible  program  structure. 
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This  is  the  generalized  equation  for  the  phase  difference  between  two 
counter-rotating  beams  in  a  Sagnac  interferometer.  This  phase  difference 
creates  an  interference  pattern  that  can  be  observed  with  a  photo-detector. 

The  output  of  the  detector  is  typically  written  as 

I  =  IQ  (1  -  cos  A<|>). 

Therefore,  by  detecting  a  change  in  the  intensity  of  the  pattern,  the  rotation 
rate  of  the  interferometer  can  be  determined. 

For  a  heterodyne  detection  scheme,  the  phase  equation  for  a  single  frequency 
can  be  written  as 

A(j>  =  luSk*2  +  ilfs.  (r^  -  n2)  +  (u^  -  mj). 

2c  X 

One  can  clearly  see  that  the  two  terms  containing  the  index  of  refraction  are 

fundamental  fiber  error  source  terms  which  exist  in  all  formulations.  The 

second  is  a  true  bias  error  generated  by  a  difference  in  the  index  of 

refraction  for  the  two  beams.  For  a  typical  gyro,  relative  stability  of  the 

12 

two  indices  must  be  better  than  5  x  10  to  reduce  this  term  to  the  level 
of  1  deg/hr  noise. 

The  natural  sensitivity  of  the  cosine  detection  function  is  maximized  for  an 
argument  of  tt/2.  The  standard  homodyne  technique  adjusts  the  (m^  -  m2) 
term  in  the  previous  equation  to  achieve  this  phase  separation  between 
counter-rotating  beams  tn  optimize  performance.  The  required  phase  shift  can 
be  accomplished  by  an  acousto-optic  modulator  which  provides  a  different 
frequency  for  the  counter-rotating  beams.  This  causes  the  counter-clockwise 
and  clockwise  path  to  differ.  An  alternative  way  to  introduce  the  x/2  phase 
shift  is  a  nonreciprocal  device  such  as  a  Faraday  cell. 


The  Faraday  cell  is  constructed  of  a  material  which  exhibits  a  strong 
magneto-optic  effect  to  gain  the  required  phase  shift  with  a  small  current. 
The  Faraday  cell  stability,  however,  is  a  direct  error  term  in  the  phase 
equation.  This  requires  that  overall  instability  of  the  cell  be  much  less 
than  the  required  minimum  detectable  rotation  rate.  To  first  order,  the 
change  of  the  light  phase  due  to  the  Faraday  effect  can  be  written  as: 

0  =  HLV 

where  H  is  an  applied  magnetic  field,  L  is  the  length  of  the  cell,  and  v  is 
the  Verdet  constant.  The  effect  in  0  can  be  written  as: 

40  1  4H  iV 
°  L  H  V 

If  the  instrument  is  required  to  measure  10-5  radians  and  the  Faraday  cell 
is  producing  a  ir/2  phase  shift,  then 

js5x  ltf6. 

This  sets  the  bounds  on  the  stability  of  V,  H,  and  L.  Typical  thermal 
coefficients  on  the  Verdet  constant  are  1  x  10“4  per  °C.  This  is  nearly  2 

_5 

orders  of  magnitude  worse  than  required  for  10  radian  stability.  The 
phase  shift  stability  of  the  Faraday  cell  represents  a  major  error  in  a 
homodyne  detection  scheme.  The  detector  intensity  for  any  homodyne  detection 
scheme,  ignoring  fiber  effects,  can  be  written  as 

I  =  IQ(t)  [1  -  cos  (kfl  +  0(t))] 

where  IQ  and  0  are  varied  in  synchronous  detection  schemes.  Clearly,  if 
any  unrecognized  change  in  I  or  0  occurs,  the  net  result  would  be  a  gyro 
bias  error.  Bias  errors  due  to  intensity  fluctuations  are  generally  removed 
by  using  two  detectors  in  portions  of  the  interference  pattern  which  are  out 
of  phase  by  7r.  The  interference  patterns  at  the  detectors  can  be  expressed 
as 
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Ix  =  I0(t)  Cl  -  cos  (k a  *  0( t ) ) ] 

I2  =  IQ(t)  [1  -  cos  (kfl  +  0(t)  +  it)]. 

Results  can  be  summed  and  differenced  to  remove  the  laser  power  term.  For 
example,  let 


The  gyro  output  is  optimized  when  0  =  tt/2.  Expanding  0  to  be  0(t)  = 

-ir/2  +6  (t),  where  6(t)  is  the  small  time  dependent  fluctuations  of  the 
phase  modulator,  the  output  Z  becomes 

Z  *  6 (t)  +  k «. 

The  fluctuations  of  the  phase  modulator  appear  directly  as  a  gyro  bias  error, 
causing  an  erroneous  rate  measurement.  The  dilemma  is  that  there  appears  to 
be  no  way  to  remove  this  bias  error  in  a  standard  homodyne  measurement 
technique. 

DUPL-POLARIZATION  CONCEPT 

The  Honeywell-designed  dual-polarization  fiber  optic  gyro  is  a  modification 
and  expansion  of  the  standard  homodyne  measurement.  The  basic  principle  is  to 
simply  operate  two  independent  rate  gyros  simultaneously  in  the  same  fiber  and 
components.  Each  of  the  gyros  is  a  complete  homodyne  rate  gyro;  thus,  the 
errors  caused  by  heterodyning  are  avoided.  If  the  phase  modulator  can  be  made 
to  act  oppositely  for  the  two  gyros,  then  it  is  possible  to  combine  the 
outputs  of  the  two  gyros  so  that  the  modulator  and  its  errors  are  removed. 

This  is  in  principle  the  same  as  a  so-called  DILAG  laser  gyro.  The  two 
Independent  gyros  can  be  made  to  coexist  if  their  polarizations  are 
orthogonal.  The  desired  phase  modulation  will  be  shown  to  be  a  natural 
outgrowth  of  the  orthogonal  polarizations.  The  results  of  proper 
implementation  of  the  concept  will  be: 


o  Optimum  n/2  phase  shift  for  gyro  sensitivity 
o  2  x  increase  in  scale  factor 
o  No  first-order  length  dependent  bias  errors 
o  Cancellation  of  gyro  bias  errors  due  to  the  modulator. 

To  aia  in  understanding  the  dual-polarization  performance  with 
polarization-preserving  fibers,  let  us  define  a  coordinate  system 
convention.  For  this  discussion,  we  will  use  a  coordinate  system  tied  to 
the  light  and  aligned  with  the  orthogonal  vertical  and  horizontal  axes  of 
the  fiber.  Figure  1  shows  examples  of  linear  polarization  states 
propagating  in  opposite  directions  in  two  fibers.  The  box  diagrams 
alongside  the  fibers  are  how  the  polarization  state  appears  to  the  observer 
as  he  is  "looking  into  the  light"  with  the  direction  of  propagation  coming 
out  of  the  page.  If  the  light  is  circular  polarized,  the  direction  of 


Figure  1.  Fiber  Polarization  Convention 
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electric  vector  rotation  in  the  fiber  is  represented  by  an  arrow  in  the 
diagram  with  clockwise  rotation  indicating  right  circular  polarization. 

This  convention  will  be  used  in  following  illustrations  to  trace  the 
polarization  state  through  the  dual  polarization  fiber  gyro.  In  some 
diagrams,  other  information,  such  as  the  orientation  of  wave  plate  fast/slow 
axes  or  accumulated  phase,  will  be  presented  with  the  diagrams. 

Figure  2a  shows  a  dual  polarization  gyro  configuration  using  isotropic  or 
birefringent  fiber.  The  source  polarization  is  aligned  with  the  fast  axis 
of  the  fiber.  This  beam  is  divided  at  the  beam  splitter,  whose  axis  of 
polarization  is  such  that  the  reflected  beam  undergoes  a  it  phase  change. 

The  reflected  beam  passes  through  the  first  X/4  plate,  becomes  RCP  (right 
circular  polarized),  passes  through  the  cell,  and  undergoes  a  phase 
retardation.  The  second  X/4  plate  returns  the  RCP  light  to  the  same 
linear  polarization  as  what  entered  the  first  X/4  plate.  This  is  a  result 
of  having  the  two  plates  rotated  90  deg  with  respect  to  each  other.  The 
linear  polarized  light  enters  the  fiber  aligned  with  the  fiber  fast-axis, 
and  propagates  clockwise  around  the  loop,  causing  a  Sagnac  phase 
retardation.  The  output  from  the  fiber  is  reflected  at  the  beam  splitter, 
accumulates  another  n  change  in  phase,  and  strikes  the  detector.  A 
perfectly  isotropic  fiber  or  one  with  polarization  compensation  can  be  used 
as  well  as  a  polarization-preserving  fiber. 

The  other  half  of  the  light  from  the  source  beam  passes  through  the  beam 
splitter  and  through  the  fiber,  picking  up  a  Sagnac  phase  advance.  In  wave 
plate  #2  it  becomes  RCP,  and  since  it  is  passing  through  the  Faraday  cell  in 
the  opposite  direction,  it  now  experiences  a  phase  advance.  Wave  plate  #1 
returns  the  RCP  light  to  the  original  linear  polarization,  which  passes 
through  the  beam  splitter  and  combines  with  the  counter-rotating  beam  to 
form  an  interference  pattern  whose  phase  difference  can  be  expressed  as 

A4>1  =  2S2  +  21. 
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Figure  2a.  Dual  Polarization  Gyro  with  S-Polarized  Input  Light 
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Now  consider  Figure  2b,  which  is  the  same  interferometer  with  the 
polarization  of  the  source  rotated  by  90  degrees.  The  reflected  beam  is  not 
phase  changed.  Wave  plate  #1  creates  LCP  light  which  is  phase  advanced  in 
the  Faraday  cell.  Wave  plate  #2  returns  the  light  to  the  initial 
polarization  state;  the  rotating  fiber  coil  generates  a  phase  retardation, 
and  the  output  is  reflected  onto  the  detector.  The  transmitted  portion  of 
the  source  light  is  phase  advanced  in  the  fiber,  and  is  changed  to  LCP  light 
by  wave  plate  #2.  Because  the  LCP  light  is  passing  through  the  cell  in  the 
opposite  direction  from  the  other  beam,  it  will  receive  a  phase 
retardation,  wave  plate  #1  converts  the  output  of  the  Faraday  cell  back  to 
the  original  polarization,  which  passes  through  the  beam  splitter  to  form  an 
interference  pattern  with  the  counter-rotating  beam.  The  phase  separation 
of  this  inteferometer  can  be  expressed  as 

A$2  =  2ft  -  21. 

The  phase  output  of  the  two  polarization  situations  can  be  shown  to  be 
A<J>^  =  2ft  +  |  21  +  8  ^ 

A4>2  =  2ft  -  21  +6<J>2 

where 

ft  =  Sagnac  phase  shift  due  to  rotation 

I  =  phase  shift  introduced  by  a  Faraday  cell 

6  <j>i  2  -  non-reciprocal  fiber  effects  for  the  fast  and 
slow  fiber  polarizations. 

The  intensity  at  the  detector  is  therefore 
Ex  =  1  -  COS  (  2ft  +  (21  +  6  <t>x) ) 

E2  «  1  -  COS  ( 2ft  -  (21  -  $  *2) .) 


The  difference  of  these  intensities  can  then  be  expressed  as 
£E  =  (  ^  +  <j>2)  +  4ft  fl  +  1/2  I2  +  1/4  (  <{>2  +  <)>2) 


Notice  that,  just  as  in  the  standard  homodyne  technique,  the  fiber 
non- reciprocal  terms  5<J>^  +  6  4»2  appear  directly  as  a  bias  error. 

However,  the  phase  modulator  instability  enters  as  a  scale  factor  error.  The 
significance  of  this  can  be  demonstrated  in  a  simple  example.  Consider  a 
homodyne  measurement  with  a  perfect  fiber  of  length  1  km  in  a  0.1m  diameter 
loop,  and  a  perfect  laser  source  at  0.8p.  The  output  phase  shift  would  be 

z  .  a  +  6 
1  Xc 

where  ft  is  the  input  rate  and  $  is  the  phase  modulator  error.  For  an  input 
rate  of  10  deg/hr  (4.8  x  10-5  rad/sec), 

=  (1.27  x  10"4  +  6)  rad. 

For  the  same  measurement  in  the  dual-polarization  mode, 


or 

Z2  =  2.54  x  10"4  (l  +  62/2  )  rad. 

Now,  if  an  undetected  modulation  shift  of  1  x  10“4  rad  were  present,  the 
output  phase  shifts  would  be 

Zx  =  2.27  x  10-4  rad 

z2  =  2.54  X  10"4  +  1.27  X  10"12  rad. 


17 


For  this  10  deg/hr  input  rate,  the  conventional  homodyne  gyro  would  read  out  a 
rate  of  17.9  deg/hr,  while  the  dual-polarization  gyro  would  read  out  10.006 
deg/hr.  The  result  clearly  shows  that  the  modulator  stability  in  the 
dual-polarization  scheme  is  no  longer  the  major  error  source  that  it  is  in  the 
standard  detection  schemes. 


SECTION  3 


SENSOR  MODELING  AND  COMPUTER  PROGRAMS 


INTRODUCTION 

In  order  to  model  gyro  performance,  Honeywell  has  developed  a  numerical 
interferometer  simulation  package.  This  capability  has  provided  a  unique 
insight  into  gyro  component  effects.  The  modeling  results  presented  in  this 
section  quantify: 

o  dual  polarization  gyro  performance 

o  fiber  depolarizer  design  issues 

o  the  effect  of  birefringent  fiber  splices 
o  significant  gyro  errors  from  distributed  3  dB  couplers. 

The  most  notable  feature  of  this  package  is  its  applicability  to  arbitrarily 
large  or  complex  systems  with  dual  polarization,  multiple  paths,  and  many 
optical  components.  In  addition,  the  optical  system  is  described  by  input 
lines  at  execution  time,  rather  than  as  a  part  of  the  program  logic.  In 
addition  to  the  conventional  optical  elements,  the  input  lines  specify  input 
excitation,  intensity  detectors  and  polarization  analyzers,  and  direct  the 
type  of  disposition  of  the  detector  output  results. 

In  the  research  and  development  of  fiber  gyros  there  is  a  need  for  a 
capability  to  analyze  and  simulate  the  optical  system.  In  order  to  predict 
performance,  sensitivity,  and  the  effect  of  changing  parameters  and  imperfect 
components,  some  form  of  quantitative  calculations  would  be  very  useful. 
Analytical  methods  quickly  become  intractable  as  the  complexity  of  the  system 
increases.  For  even  the  simplest  fiber  gyro  configuration,  for  several 
optical  paths,  and  for  a  few  optical  components,  it  is  practically  impossible 


to  reduce  a  complete  analytical  description  to  a  useful  form.  An  accurate 
numerical  simulation  would  be  invaluable  for  low  cost  exploration  of  different 
designs,  sensitivity  calculation,  and  performance  prediction. 

Honeywell  has  been  developing  a  program  for  such  a  simulation,  and  we  believe 
that  it  will  be  a  useful  tool  in  such  work.  The  novel  features  of  the  program 
are: 

1.  It  is  applicable  to  arbitrarily  large  systems,  with  multiple 
optical  paths  and  a  large  number  of  optical  components. 

'  2.  The  optical  system  is  described  by  a  sequence  of  input  lines, 

not  by  the  program  logic.  The  program  does  not  need  to  be 
changed  for  each  optical  structure  to  be  studied. 

3.  Different  types  of  "experiments"  can  be  run,  each  one  again 
determined  by  single  input  lines. 

4.  The  program  is  easily  expandable  to  include  additional  optical 
elements  or  experiments. 

METHOD:  JONES  VECTOR-MATRIX  REPRESENTATION 

The  optical  system  is  modeled  using  the  well-known  Jones  vector-matrix 
representation  for  the  wave  amplitude  propagation  along  multiple  optical 
paths.  The  effect  of  optical  components  on  the  light  is  represented  by  matrix 
operations.  Wave  amplitudes  are  represented  by  complex  two-vectors  for  the 
horizontal  and  vertical  (x  and  y)  components,  here  denoted  as  Ax  and  Ay. 

Optical  components  are  represented  by  2  x  2  complex  matrices  whose  diagonal 
elements  may  be  viewed  as  propagation  factors  for  each  of  the  x  and  y  wave 
components;  the  off-diagonal  elements  are  coupling  or  conversion  coefficients 
between  the  Ax  and  Ay.  Standard  matrices  are  well-known  for  most  conventional 
optical  components  such  as  polarizers,  isotropic  or  birefringent  path  lengths, 
wave  plates,  retarders,  rotators,  etc. 
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For  modeling  interferometric  systems  it  is  essential  to  retain,  along  each 
path,  information  on  the  total  wave  amplitude  (power)  and  the  absolute  phase 
along  the  path.  This  is  done  by  including  a  common  complex  propagation 
factor  associated  with  each  two-vector.  Conventionally,  the  Jones 
representation  assumes  normalized  vectors  and  matrices  such  that  the  total 
power  in  a  path  is  irrelevant,  and  any  phase  factor  common  to  both  Ax  and  Ay 
is  assumed  to  be  factored  out.  With  such  an  assumption  it  is  not  possible 
to  compare  (add,  combine,  or  couple)  the  waves  in  two  paths.  Thus,  the 
present  method  may  be  described  as  the  Jones  vector-matrix  method,  enhanced 
by  preservation  of  the  absolute  phase  and  total  amplitude  along  each  path. 

Beam  splitters,  combiners,  and  distributed  directional  couplers  are  also 
necessary  optical  components.  These  are  slightly  less  familiar  to  the  usual 
descriptions  of  the  Jones  method,  but  the  form  of  the  required  2x2  and 
4x4  matrices  is  well-known. 

The  vector-matrix  representation  is  ideally  suited  to  numerical  work  and  to 
computer  programs.  The  mathematical  quantities  are  easily  calculated  for 
any  component  and  values  of  the  parameters  which  describe  it,  e.g.,  path 
lengths,  orientation,  reflection  and  transmission  coefficients,  coupling, 
etc.  The  complexity  of  a  system  lies  in  the  large  number  of  its  components 
and  the  multiplicity  of  interconnecting  paths,  rather  than  in  the  basic 
components  themselves.  This  problem  can  be  handled  by  a  good  computer 
program  structure. 

PROGRAM  STRUCTURE 


1 

I 

I 
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Logically  there  are  three  main  parts  to  the  program.  The  first  part  reads 
the  input  lines,  one  per  optical  component,  and  builds  up  to  a  description 
of  the  optical  system.  Pointer  vectors  and  lists  are  built  up  which 
describe  the  operations  to  be  carried-out,  in  which  order,  and  on  which 
paths.  The  operations  are  simply  the  multiplication  of  the  complex  wave 
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amplitude  of  each  path  by  the  2x2  matrix.  This  phase  of  the  program 
execution  may  be  described  as  a  "compilation"  of  the  optical  system,  to  be 
executed  later  any  number  of  times. 

The  compilation  of  the  system  description  from  the  input  lines  is  done  only 
once.  But  execution  for  numerical  output  may  be  carried  out  any  nunber  of 
times  to  give  the  output  data  points  vs  some  changing  parameter  such  as 
wavelength,  a  rotation  rate,  or  path  length,  etc. 

The  heart  of  the  program  is  a  subroutine  containing  a  large  number  of  small 
blocks  of  code  for  matrices  and  multiplications.  A  block  of  lines,  from  a 
few  to  perhaps  20,  exists  for  each  possible  component.  During  execution, 
for  any  fixed  set  of  parameter  values,  and  progressing  from  input  to  output, 
the  logical  flow  is  to  skip  to  each  block  in  order,  executing  only  those 
operations  called  for.  This  is  done  very  efficiently  following  the  pointer 
lists  made  up  during  the  build  mode.  It  is  done  repeatedly  with  changing 
parameter  values  during  any  one  "experiment." 

The  third  part  of  the  program  is  another  subroutine  which  conducts  the 
"experiments"  or  simulation  runs.  It  reads  additional  input  lines,  one  or 
more  per  experiment,  which  describe  the  parameters  to  be  varied  and  the 
values  to  be  used,  and  then  calculates  the  "data  points."  The  detector 
outputs  are  printed  or  stored  for  each  data  point.  Any  number  of  different 
experiment  runs  can  be  carried  out  by  successive  input  lines,  all  for  the 
same  optical  system,  without  repeating  the  system  description. 

PROGRAM  INPUT:  AVAILABLE  COMPONENTS  AND  EXPERIMENTS 

The  structure  of  the  optical  system  and  the  types  of  numerical  experiments 
to  be  simulated  are  controlled  by  a  list  of  input  lines.  One  type  of  input 
line  defines  an  optical  component.  It  specifies:  a  mnemonic  for  the 
component;  the  path  or  paths  of  which  it  is  a  part;  and  up  to  four  real 
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parameters  of  the  component.  A  different  type  of  input  line  chooses  the 
type  of  numerical  experiment  to  be  performed:  it  specifies  the  particular 
optical  component(s)  and  parameter(s)  to  be  varied,  and  defines  the  list  or 
range  of  values  to  be  used,  as  well  as  the  disposition  of  the  generated 
output. 

Standard  optical  components  include:  path  lengths,  both  isotropic  and 
birefringent  with  arbitrary  orientation;  polarizers  for  linear,  circular, 
and  arbitrary  elliptical  polarization;  Faraday  cell;  beam  splitters  and 
combiners,  both  polarizing  and  nonpolarizing;  polarization  rotators; 
waveplates  and  retarders  of  quarter-,  half-  and  arbitrary  length  and  any 
orientation;  distributed  couplers  in  terms  of  coupling  coefficient, 
coupling  length  and  phase  mismatch;  and  beam  splitters'  combiners. 

The  list  of  available  optical  components  includes  non-standard  components 
such  as  detectors  of  total  intensity  or  intensity  in  specific  polarizations, 
and  a  Stokes  parameter  detector,  allowing  the  polarization  state  to  be 
plotted  on  the  Poincare  sphere.  These  detectors,  which  may  be  positioned  at 
any  point  in  the  optical  system,  do  not  affect  beam  propagation.  Other 
non-standard  components  are  a  non-reciprocal  (rotation)  simulator  and  a 
"component"  which  specifies  the  excitation  state  of  the  input  beam. 

Additional  pseudo-components  are  available  which  can  duplicate  or 
symmetrically  reflect  a  previously  input  component  (or  set  of  components)  to 
a  position  later  in  the  chain.  One  novel  input  "component"  with  mnemonic 
FLOP  can  unfold  a  gyro  structure.  Any  number  of  components  preceding  the 
midpoint  of  the  fiber-coil  or  circular  path  are  reflected  about  the 
midpoint,  with  path  numbers  interchanged  and  sense  of  propagation  reversed. 
This  economizes  the  input  required,  and  guarantees  duplication  of  the 
components  and  their  parameters  for  the  counterpropagation  paths.  The 
latter  guarantee  is  important  in  experiments  which  vary  a  parameter  of  a 
component  in  the  counter-propagating  paths  of  the  system. 
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Various  types  of  experiments  may  be  performed  by  varying:  1)  any  one 
parameter  of  one  or  several  duplicated  components;  2)  the  relative 
frequency  of  the  input  light  beam;  and  3)  any  two  parameters  in  one  (or 
duplicated)  components.  Parameters  may  be  varied  by  specifying  a  list  of 
values,  or  by  incrementing  over  a  range. 

EXAMPLES 

These  examples  are  presented  in  this  section  following  the  Discussion  to 
illustrate  model  capabilities: 

1.  Dual  Polarization  Gyro. 

2.  Narrow-band  Birefringent  Depolarizer. 

3.  Four  Sections  of  Birefringent  Fiber. 

4.  Basic  All-Fiber  Gyro,  Effect  of  Coupler  Drift. 

These  results  are  important  because  they  represent  a  detailed  examination 
of  fundamental  gyro  error  sources.  The  first  example  is  central  to  this 
contract  program;  the  others  are  important  to  recently-developed 
high-performance  gyro  configurations. 

DISCUSSION 

The  program  is  just  reaching  fruition,  so  these  examples  only  hint  at  the 
principal  power  of  the  program — the  ability  to  do  simulations  of 
arbitrarily  complicated  optical  systems  with  multiple  paths  and  many 
optical  components. 


The  limitations  of  the  program  are  those  imposed  by  the  Jones  method. 

Only  monochromatic  steady  state  can  be  treated  at  present.  Consequently, 
the  program  is  not  usable  for  time-modulated  systems  which  are  becoming 
important  for  fiber  gyros.  For  quasi-monochromatic  excitation  (e.g., 
semiconductor  lasers,  LEDs,  etc.),  it  is,  in  fact,  possible  to  extend  the 
Jones  formulation  and  the  program.  It  is  only  necessary  to  replace  the 
complex  two-vector  of  the  wave  amplitudes  by  their  coherency  matrix. 

This  is  a  trivial  change  in  the  program;  the  same  Jones  matrices  are 
used  for  the  optical  components.  It  is  conceivable  that  the  program  can 
be  extended  to  treat  the  Lime-modulated  cases,  but  this  would  probably 
require  considerably  more  complexity  to  keep  track  of  transit  time 
effects.  The  method  is  also  not  applicable  to  resonator  structures; 
those  become  boundary  value  problems. 

In  actual  practice,  an  important  limitation  will  be  lack  of  adequate 
quantitative  data  to  use  with  the  program  for  realistic  modeling.  Many 
optical  components  are  less  than  ideal,  and  too  little  is  known  about  the 
parameters.  Examples  are  the  polarization  effects  in  beam  splitters  and 
distributed  couplers,  and  twist  or  misalignment  in  birefringent  fibers. 
The  modeling  can  be  only  as  good  as  the  optical  parameters  are  well-known. 

Example  1.  Dual  Polarization  Gyro 

Figure  3  shows  the  gyro  structures  as  it  can  be  modeled,  both  in  its 
folded  and  unfolded  form.  All  significant  optical  components  are 
included.  Input  excitation  is  at  +45  deg  polarization,  and  separate 
detectors  for  x  and  y  polarizations  are  used  at  the  output.  The 
wavelength  used  was  0.633.10”3  mm. 
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Figure  3.  Dual  Polarization  Gyro  as  Modeled.  Top,  physical  form; 
bottom,  unfolded  form  simulated  by  computer  model. 


Polarization-preserving  fiber  with  a  An  =  0.15  x  10"^  was  used  for 
the  coil  of  100  meters.  Because  a  model  is  not  yet  included  for  a  path 
which  is  both  birefringent  and  rotation  sensitive,  the  length  was 
partitioned  into  alternating  segments  of  birefringent  fiber  and 
nonreciprocal  fiber,  for  a  total  of  50  meters  of  each  type  (with  twice 
the  retardation  coefficients  appropriate  for  a  100m  length).  To  simulate 
fiber  twist  or  misalignments,  the  orientation  of  the  3  birefringent 
sections  could  be  set  at  +  a  few  degrees. 

Figure  4a  shows  the  output  from  the  x  and  y  polarization  detectors  vs 
rotation  rate  in  deg/sec.  Scale  factor  and  drive  for  the  Faraday  cell 
were  adjusted  for  the  quadrature  operating  point,  and  the  fast  axis  of 
all  the  birefringent  fiber  sections  was  aligned  to  0  deg.  The  output  is 
as  expected,  with  the  rate  signal  from  the  two  detectors  having  opposite 
phase  with  respect  to  the  direction  of  orientation.  (The  two 
polarizations  are  biased  to  complementary  quadrature  points.)  An 
arbitrary  scale  factor  was  used  to  make  this  model  sensitivity  agree  with 
the  experimentally-observed  80  deg/sec  for  2tt  phase  shift  in  the 
detector  outputs. 

To  simulate  possible  angular  misalignment  of  the  fiber,  the  three 
sections  were  then  oriented  at  +5,  +7,  and  -8  deg  respectively,  with  the 
results  shown  in  Figure  4b.  Most  notable  is  the  shift  from  equal 
'detector  output,  appearing  as  a  drift  to  an  equivalent  +10  deg/sec 
rotation  rate.  Also,  the  curves  are  no  longer  symmetric — most  notably, 
the  curve  for  the  x  polarized  output. 
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A  different  experiment  is  shown  in  Figure  5,  where  the  scale  factor  of 
the  Faraday  cell  (Verdet  constant)  was  varied  over  a  relative  range  of 
+5  *  10“\  This  was  intended  to  simulate  the  effects  of  a  temperature 
variation  of,  say,  +5  deg,  assuming  a  temperature  coefficient  of  10_4/deg. 
This  demonstrates  the  effects  of  small  changes,  as  evidenced  by  the  expanded 
scales  needed.  The  rotation  rate  was  0.002  deg/sec,  and  bias  was  at 
quadrature  at  the  0.15  point  on  the  horizontal  axis,  resulting  in  0.5  output 
for  both  detectors.  The  change  in  detector  output  with  this  "drift"  in  the 
Faraday  cell  was  identical;  the  two  curves  have  a  constant  separation. 

With  respect  to  the  sense  of  direction  of  rotation,  or  an  apparent  rotation, 
the  sense  of  detector  drift  is  opposite  for  the  two  outputs.  The  rotation 
signal,  the  difference  between  the  two  detectors,  does  not  drift.  This 
directly  demonstrates  the  advantage  of  the  dual  polarized  gyro. 

Figures  6a  and  6b  show  similar  results  in  the  presence  of  the  birefringence 
misalignment  experiment  previously  shown  in  Figure  4b.  We  used  the  same 
change  and  horizontal  scale  for  the  Faraday  cell  scale  factor  as  in  Figure 
5.  Note  the  zero  rotation  rate  points  in  Figure  4b:  at  about  0.19  for  the 
x  polarization,  and  about  0.84  for  the  y  polarization.  These  points  vs  the 
drift  in  the  bias  point  are  shown  respectively  in  Figures  6a  and  6b.  Note 
the  expanded  scales  on  the  left  for  detector  output.  The  slopes  of  the  two 
curves  are  in  the  same  direction,  but  because  the  two  are  on  such  different 
parts  of  their  respective  curves,  the  slopes  are  very  different.  The 
important  point  here  is  that  the  misorientation  at  several  "joints"  in  the 
highly  birefringent  fiber  has  introducted  mixing  of  the  x  and  y  polarized 
energy.  The  consequences  are  quite  unpredictable. 
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Figure  6a.  Dual  Polarization  Output  vs  Faraday 
Cell  Shift  for  Y  Polarizations 
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Figure  6b.  Dual  Polarization  Output  vs  Faraday 
Cell  Shift  for  X  Polarizations 
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These  are  preliminary  results  and  are  not  fully  understood.  The  degree  of 
misalignment  is  perhaps  rather  large,  and  the  abrupt  change  is  not  really 
equivalent  to  fiber  twist. 

Example  2.  Narrow-Band  Birefrlngent  Depolarizer 

Rather  than  preserving  a  polarization  state,  it  is  sometimes  desirable  to 
depolarize  a  light  source.  This  has  meaning  only  for  non-monochromatic 
light,  and  the  bandwidth  over  which  depolarization  can  be  obtained  is  an 
important  question.  As  an  example,  in  the  use  of  superluminescent  diodes  as 
sources  for  fiber  gyros,  a  depolarizer  over  a  1%  bandwidth  is  desirable. 

Figure  7  shows  a  possible  depoarizer  using  two  lengths  of  quartz.  For 

narrow-band  performance,  the  pieces  need  to  be  rather  long,  depending  on  the 

degree  of  birefringence.  One  example  given  in  the  literature  states  that 

o 

the  longer  piece  should  be  twice  as  long  as  the  short  piece.  It  can  be 
shown  that  this  is  a  poor  choice;  rather,  the  two  lengths  should  have  a 
ratio  not  near  the  ratio  of  small  integers.  The  total  length,  in  number  of 
retardation  wavelengths  at  center  band,  should  be  many  times  the  inverse 
relative  bandwidth. 

The  Poincare  sphere  is  very  useful  to  exhibit  the  spreading  of  the 
polarization  states  by  a  depolarizer.3  In  a  good  depolarizer,  as  the 
frequency  of  the  light  is  varied  within  the  bandwidth,  the  state  of 
polarization  is  redistributed  widely  and  uniformly  over  the  surface  of  the 
sphere.  The  lower  part  of  Figure  7  shows  the  points  on  the  sphere  which 
were  obtained  using  the  computer  program.  (This  model  requires  a  simple 
list  of  about  10  input  lines.)  The  relative  frequency  was  varied  over  the 
range  of  0.995  to  1.005  using  100  points.  That  is,  the  total  spectral  range 
displayed  is  one  percent,  and  the  individual  points  shown  are  spaced  at 
10 "4.  Good  distribution  over  the  sphere  is  evident,  although  there  is  a 
somewhat  greater  number  of  points  (loops)  near  the  +45  deg  linear  polarized 
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Figure  7.  Narrow-Band  Two-State  Birefringent  Depolarizer. 
Locus  of  polarization  states  on  Poincare 
sphere  over  a  1*  spectral  range. 


points  P  and  Q.  Different  choices  of  parameters  for  best  distribution  on 
the  sphere  are  very  easily  explored  using  the  computer  program.  Design 
iterations  are  efficient  because  of  the  flexible  program  structure. 


Example  3.  Four  Sections  of  Blrefringent  Fiber 

The  results  from  a  simple  example  of  birefringent  fiber  sections  misoriented 
in  angle  at  their  junctions  are  shown  in  Figures  8  and  9.  Four  sections  of 
fiber  with  a  linear  birefringence  of  6  n  =  2  •  10~7  were  used  (beat 
length  of  2.5  meter  at  X  =  0.5  •  10~3mm).  Lengths  were  9.5,  9.09, 

9.72,  and  9.45  meters;  and  the  orientations  (fast  axis)  were  all  nominally 
45  deg,  but  actually  misoriented  to  41,  47,  50,  and  45  degs  respectively. 
Input  light  polarized  at  45  deg.  is  nominally  aligned  with  the  fast  axis  of 
the  fiber,  and  should  emerge  as  45  deg  polarized  light.  But  the  fiber 
misorientation  will  result  in  deviations  from  pure  45  deg.  polarization,  as 
is  clearly  shown  in  the  two  figures  where  the  relative  frequency  was  varied 
over  the  range  0.95  (0.001). 05. 

Figure  8  shows  the  output  from  x  and  y  polarized  detectors.  The  sum  of  the 
two  is  constant  (energy  conservation);  and  where  the  two  are  equal,  the 
azimuth  of  the  polarization  is  45  deg,  though  there  will  generally  be  some 
ellipticity.  More  clearly  displaying  the  departure  from  the  +45  deg. 
polarization  state  is  the  locus  shown  on  the  Poincare  sphere  in  Figure  9. 

The  locus  spirals  around  the  point  P,  or  Stokes  parameter  axis  S2.  Over  a 
wider  frequency  band  the  spiral  path  repeats  periodically.  Projecting  this 
locus  onto  the  SI  axis  from  H  to  V  gives  the  resolution  into  x  and  y 
polarization  shown  in  the  previous  figure. 


35 


Example  4.  Temperature-Induced  Bias  Errors  in  Gyros  Using  Distributed 
Couplers 

The  basic  all-fiber  gyro  consists  of  a  distributed  3  dB  fiber  coupler  used 
for  a  beam  splitter/combiner  plus  the  rotation  sensing  fiber  coil.  This  may 
be  modeled  quite  simply,  with  a  single  non-reciprocal  phase  retarder  in  the 
path  to  establish  the  quadrature  operating  point.  Two  unfolded  paths  are 
used  for  the  clockwise  and  counterclockwise  beams  respectively,  with  a 
single  detector  for  the  output  of  each  path,  giving  complementary  rotation 
signals.  Only  a  single  polarization  is  used. 

An  experiment  of  interest  is  how  the  output  signal  might  drift  with,  say, 
temperature  changes  in  the  distributed  coupler.  It  is  not  obvious  whether 
the  drift  would  be  equivalent  to,  or  be  interpreted  as,  a  rotation  signal, 
but  this  is  easily  modeled.  Figure  10  shows  the  results  of  varying  the 
coupling  coefficient  of  the  distributed  coupler  (this  is  equivalent  to 
varying  the  effective  coupler  length  or  the  frequency  of  the  light).  The 
consequence  is  that  the  light  is  not  equally  divided  between  the  two  paths. 
It  is  the  coupling  coefficient  which  is  varied  here,  not  the  power  division 
ratio.  Because  of  the  nature  of  the  distributed  coupler,  a  phase  shift  also 
accompanies  the  power  division  change;  this  is  a  phase  change  that  appears 
as  a  gyro  bias  error. 

With  the  choices  of  wavelength  and  coupler  length  used,  equal  power  division 
required  a  coefficient  of  0.25  x  10“5  per  unit  length.  This  point  is 
shown  circled  in  Figure  10.  At  quadrature  and  zero  rate,  the  detectors  for 
the  two  outputs  are  equal.  Along  the  horizontal  axis,  the  coupling 
coefficient  is  varied  over  a  wide  range.  The  detector  outputs  vary  in  a 
manner  indistinguishable  from  a  rotation  rate.  To  first  order,  very  near 
the  tt/2  point,  the  effect  vanishes;  but  the  second  order  effect,  for 
changes  of  even  one  percent  in  the  coefficient,  still  can  be  large  compared 
to  the  rotation  rates  to  be  measured. 


38 


For  typical  gyro  parameters  and  a  few  degrees  temperature  change,  this  error 
term  amounts  to  hundreds  of  degrees  per  hour.  The  magnitude  of  this  effect 
has  not  generally  been  recognized.  It  may  present  a  fundamental  limit  to 
fiber  optic  gyro  performance. 


SECTION  4 


EXPERIMENTAL  RESULTS 


A  significant  portion  of  the  objective  of  this  contract  was  to  determine  the 
operating  parameters  of  the  dual  polarization  gyro.  The  most  fundamental 
part  of  this  is  the  definition  of  concept  feasibility.  This  includes 
demonstration  of  the  experiment  with  measurements  to  evaluate  performance  of 
the  gyro  at  quadrature,  and  especially  to  observe  the  effect  of  rotation  and 
Faraday  cell  bias  on  the  two  polarizations. 

The  most  important  result,  described  later  in  this  section,  was 
demonstration  of  dual  polarization  gyro  performance  with  a  changing  phase 
bias  element  environment.  Gyro  bias  is  shown  to  be  an  order  of  magnitude 
more  stable  in  the  dual  polarization  configuration  than  with  only  one 
polarization  operating. 

GYRO  COMPONENTS 

The  experimental  set-up  used  is  shown  in  Figure  11.  Here  optical  components 
are  arranged  on  a  plate  which  has  been  mounted  on  a  rate  table.  Fiber  is 
wound  around  the  periphery  of  the  plate.  In  this  experiment  approximately 
150  meters  of  fiber  are  used  on  a  0.9  meter  diameter  circle.  A  Tropel  100 
helium-neon  light  source  is  used  for  the  initial  experiment.  This  laser 
produced  noise  due  to  light  feedback  from  the  experiment.  In  later 
arrangements  a  Faraday  cell  isolator  was  used  to  reduce  the  returning  light, 
thereby  reducing  oscillations  in  the  laser.  Laser  output  power  is  sampled, 
as  shown,  by  detector  3.  This  power  measurement  and  others  at  various 
places  on  the  table  are  connected  to  a  dedicated  computer  facility  used  to 
record  data,  analyze  performance,  and  plot  gyro  characteristics. 


Figure  11.  Gyro  Component  Arrangement  or  Rate  Table 


The  50%  beamsplitter  which  allows  coupling  light  to  both  fiber  ends  is  made 
from  a  quartz  substrate  and  designed  to  work  at  near  normal  incidence, 
minimizing  polarization-dependent  effects.  A  second  Faraday  cell  with 
associated  quarter  wave  plates  provides  the  non-reciprocal  element  needed 
for  dual  polarization  gyro  biasing.  Quarter  wave  plates  on  either  end  of 
this  Faraday  cell  are  aligned  with  a  90°  difference  with  respect  to  each 
other. 

Fiber  end  couplers  consist  of  a  low  strain  microscope  objective  and  a  3-axis 
differential  screw  translation  stage.  One  end  of  the  fiber  passes  through  a 
polarizer4  and,  for  some  configurations,  it  has  a  fiber  stretcher  for 
decorrelation  of  fiber  scattering.  Because  the  gyro  runs  two  polarizations 
simultaneously,  it  is  necessary  to  have  detectors  arranged  for  each;-  in 
this  case,  a  polarizing  beam  splitter  is  used  for  signal  separation  at  the 
output  of  the  gyro. 

The  optical  components  on  the  rate  table  are  connected  through  a  slip  ring 
assembly  to  a  dedicated  computer  data  retrieval  system.  This  is  centered 
around  a  Fluke  model  1720  A  instrument  controller.  This  system  includes  an 
analog  scanner  which  feeds  a  remotely  controlled  6-1/2  digit  voltmeter. 

Data  is  stored  in  the  1720  A  memory,  processed,  and  then  stored  on  disk 
file.  This  portion  of  the  system  is  illustrated  in  Figure  12. 

A  function  generator  and  power  amplifier  are  also  connected  to  the  slip  ring 
assembly.  They  are  used  for  experiments  requiring  PZT  stretchers  for 
decorrelation  of  Rayleigh  scattering  within  the  fiber. 

SOFTWARE 

The  experiments  run  on  the  fiber  optic  gyro  program  used  a  Fluke  1720A  data 
controller  as  the  main  data  storage  and  analysis  device.  The  1720A  is  a  16 
bit  micro  computer  with  a  floppy  disc  storage  media.  The  language  used  is 


43 


enhanced  Basic.  All  the  special  drivers  to  control  the  IEEE-488  data  bus 
which  connects  the  outside  world  to  the  CPU  are  resident  in  the  Basic 
operating  system. 


The  routines  used  to  acquire,  analyze  and  plot  the  experimental  data  were 
generated  by  Honeywell  personnel.  They  make  full  use  of  the  1720A's 
capabilities  as  a  simple  dedicated  experiment  controller.  The  programs 
developed  are  DATRUN,  a  data  acquisition  routine;  DPDEV,  a  data  acquisition 
and  real  time  analysis  routine;  DPSTAT,  a  statistical  analysis  routine; 
RATMAX,  a  rate  data  pre-processing  routine;  and  HPL0T1,  a  data  plotting 
routine.  Details  of  these  software  packages  are  contained  in  Appendix  A. 

STATISTICS  AND  NORMALIZATION 

The  output  of  the  dual  polarization  gyro  detectors  is  a  voltage.  The 
functional  form  of  this  voltage  is 

V  =  V°  sin  Kft 

where  V°  is  the  intensity  scale  factor  of  the  output  and  K  is  the  rate  scale 
factor.  K  is  determined  uniquely  by  the  gyro  parameters  and  is  written  as 

k  =  *-LR- 

Xc 

where 

L  =  fiber  length  =  150m 

R  =  loop  radius  =  ,5m 

X  =  optical  wavelength  =  .633xl0“6m 

Q 

c  =  speed  of  light  =  3x10  m/sec. 

For  this  gyro,  KD  =  5  radians/rad/sec. 

The  intensity  scale  factor  V°  is  determined  by  the  laser  power,  fiber 
losses,  fiber  coupling  coefficients,  general  set-up,  and  detector  scale 
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factor.  Because  of  the  flexibility  of  our  experimental  set-up,  this  Vo  has 
had  the  opportunity  to  vary.  We  made  no  real  effort  to  fix  this  value 
throughout  our  tests.  We  only  required  that  it  was  known  for  the 
configuration  being  tested.  The  value  is  easily  determined  experimentally  by 
a  rate  variation  experiment  where  the  input  rate  is  varied  such  that  xn 
traverses  an  angle  greater  than  2tt. 

For  the  most  part,  Vo  =6  volts.  Some  data  plots  have  a  different  Vo.  In 
those  cases,  Vo  is  recorded  on  the  plot. 

In  order  to  determine  the  gyro  performance  in  real  units,  i.e.,  deg/hr,  we 
would  normally  have  to  invert  the  equation  V  =  Vo  sin  Kft  to  the  form 


Practically  speaking,  this  is  a  difficult  process  to  implement  and  is  not 
really  necessary  for  experiments.  The  sensitivity  of  the  gyro  can  be 
described  in  terms  of  the  standard  deviation  of  the  rate  data.  This  can  be 
written  in  terms  of  as 

°b  ‘Y/ovL5™  °v 

or,  for  a  zero  rate  drift  run, 


Because  the  dual  polarization  can  operate  at  quadrature  at  zero  rate,,  the 
voltage-to-rate  conversion  can  be  approximated  at  low  rates  by 


Therefore,  the  zero  rate  drift  test  results,  which  are  presented  in  volts,  can 
be  easily  converted  into  true  rate  values. 

EXPERIMENTAL  RESULTS 

The  dual  polarization  gyro  principle  described  in  Section  2  of  this  report  was 
verified  and  is  illustrated  in  the  following  figures.  Using  the  breadboard 
set-up  of  Figure  11,  we  were  able  to  change  optical  components  and  vary 
experimental  conditions  to  optimize  the  gyro  operation.  Figure  13  illustrates 
a  dual  polarization  gyro  output  showing  the  horizontal  and  vertical  detector 
signals.  This  plot  shows  several  of  the  unique  features  of  dual  polarization 
configurations. 

At  zero-rate  the  two  outputs  are  observed  to  be  at  the  quadrature  point 
(V2),  and  of  opposite  slope.  This  situation  is  obtained  with  the  gyro 
aligned  for  maximum  output.  These  two  polarization  outputs  have  been  combined 
algebraically  in  Figure  14.  Comparing  Figures  13  and  14,  one  can  note  that 
the  scale  factor  for  the  gyro  has  been  increased  by  a  factor  of  two  and  the 
noise  decreased.  The  nature  of  the  noise  is  obvious  from  re-examining  Figure 
13,  in  which  the  correlated  peaks  can  be  seen. 

NOISE  REDUCTION  TECHNIQUES 

In  order  to  decrease  the  noise  in  the  gyro  output,  we  analyzed  noise  sources. 
The  most  fruitful  portion  of  this  centered  on  laser-generated  noise.  Figure 
15  illustrates  a  typical  gyro  output  power  spectral  density  for  a  drift  run 
(constant  rate  input).  Here  the  noise  peak  at  approximately  4.0  KHz  is 
similar  to  the  power  spectral  density  noise  peak  for  the  laser  output.  The 
laser  output  power  spectral  density  is  shown  in  Figure  16.  This  information 
led  us  to  conclude  that  some  form  of  laser  isolation  was  required. 
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Figure  13.  Gyro  Output  vs.  Rotation  Rate  for  Horizontal 
and  Vertical  Gyro  Polarizations 
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Figure  14.  The  Combination  of  Two  Polarization  Outputs 
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Figure  15.  Power  Spectral  Density  of  Gyro  Output 


We  constructed  a  Faraday  rotator/polarizer  using  materials  suggested  by  the 
Naval  Research  Laboratory.5  This  isolator,  shown  in  Figure  17,  was 
instrumental  in  reducing  overall  gyro  noise  levels.  Other  techniques,  such 
as  improved  optical  component  quality  and  spatial  filtering  of  the  gyro 
input,  were  also  examined. 

These  improvements  resulted  in  a  significant  decrease  in  system  noise,  which 
allowed  backscattering  from  the  fiber  to  become  the  most  significant  noise 
source.  The  fiber  stretcher  was  implemented  to  decorrelate  backscattering, 
allowing  further  improvement  in  the  gyro  performance.  Both  frequency  swept 
(chirped)  and  continuous  sine  wave  exertation  for  the  PZT  stretcher  were 
evaluated.  The  most  effective  decorrelation  was  found  to  be  in  the  4.0  - 
4.5  KHz  region.  A  1-centimeter  diameter  PZT  cylinder  with  about  15V  was 
used  with  10  fiber  turns.  Figure  18  illustrates  the  before/after  results  in 
a  drift  run  using  the  dual  polarization  gyro  output.  The  RMS  noise  in  the 
runs  in  Figure  18  show  a  factor  of  31  improvement  with  the  decorrelator  on. 

The  gyro  output  characteristics  using  the  noise  reduction  techniques 
described  above  are  shown  in  Figure  19.  This  output  compares  favorably  with 
the  previous  output,  Figure  14.  The  scale  factor  change  in  the  two  figures 
is  a  result  of  neutral  density  filters  in  the  gyro  input  portion  of  the 
experiment.  At  the  time  these  data  were  taken,  the  input  power  was 
approximately  100  microwatts. 

A  long-term  drift  run  is  shown  in  Figure  20.  Here  the  data  are  plotted  with 
an  expanded  vertical  scale,  allowing  more  detailed  examination  of  the 
noise.  Aside  from  a  turn-on  transient,  the  data  are  relatively  stable.  The 
drift  indicated  by  the  portion  of  the  plot  after  the  stabilization  is 
approximately  50  deg/hr/hr.  This  is  a  surprisingly  good  bias  stability 
value  for  such  a  flexible  experimental  set-up. 


52 


Figure  17.  Faraday  Rotator 
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Figure  19.  Dual  Polarization  Gyro  Output  After  Noise  Reduction 
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Figure  20.  Long-Term  Drift  Run  with  Expanded  Vertical  Scale 
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The  most  important  objective  of  this  experimental  effort  was  to  demonstrate 
dual  polarization  gyro  independence  from  phase  modulator  errors.  It  was 
hoped  that  the  analytical  predictions  of  Section  2  and  the  numerical 
modeling  of  Section  3  could  be  verified  experimentally.  To  accomplish  this 
it  would  be  necessary  to  measure  gyro  bias  while  changing  phase  modulator 
settings.  The  principal  factor  in  phase  modulator  stability  is 
environmental  temperature;  this  parameter  was  varied  over  a  time  period  of 
about  one  hour,  and  the  gyro  drift  recorded. 

Figure  21  shows  the  temperature  change  of  the  Faraday  cell  housing  as  a 
function  of  time.  In  this  example  it  was  heated  to  about  33°C  and  allowed 
to  cool.  During  the  experiment,  laser  output  power  was  also  observed  to 
fluctuate.  This  is  plotted  in  the  lower  curve  of  Figure  21.  A  conventional 
gyro  would  be  expected  to  respond  (drift)  in  response  to  these 
perturbations.  This  is,  in  fact,  what  is  observed  in  Figure  22.  The  curves 
labeled  0  and  2  are  the  two  single  polarization  outputs  of  a  gyro  plotted  on 
the  same  time  scale  shown  in  the  previous  figure. 

Curve  1  in  Figure  22  shows  the  dual  polarization  performance  of  this  gyro. 

It  is  an  order  of  magnitude  more  stable  than  the  conventional  (single) 
polarization  gyro.  This  confirms  the  expectation  of  better  drift 
performance  for  the  dual  polarization  gyro. 


POWER  DRIFT) 


Figure  21.  Faraday  Cell  Temperature  and  Laser  Output  Power 
Fluctuations  vs  Time  During  Drift  Measurements 
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Figure  22.  Single  Polarization  (curves  0,  2)  and  Dual 

Polarization  (curve  1)  Output  Showing  an  Order 
of  Magnitude  Stability  Improvement  for  Dual 
Polarization. 
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APPENDIX  A 


SOFTWARE  DETAILS 


Programs  used  for  Fluke  1720A  programming  are  included  in  this  report  to 
allow  their  use  by  other  researchers,  and  to  quantify  the  statistical 
measures  of  performance  used  in  the  experimental  portions  of  this  report. 
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DATRUN 

This  routine  is  a  general-purpose,  fast  data  acquisition  routine  that  reads 
and  stores  up  to  10  channels  of  6-1/2  digit  analog  data  at  a  rate  up  to  200 
channels  per  second.  The  routine  sets  up  the  7600-scanner/DVM-subsystem  in 
its  fastest  acquisition  mode.  The  number  of  channels  read  and  the  speed  at 
which  they  will  be  read  are  under  user  control  at  run  time.  On  command  the 
program  acquires  data  and  stores  them  in  temporary  storage  until  the 
operator  terminates  the  test  or  the  storage  space  is  exhausted.  This  read 
and  store  operation  allows  for  the  high-speed  operation.  At  completion  of 
the  test,  stored  data  are  converted  to  real  numbers  and  stored  on  the 
diskette  in  a  permanent  file.  The  data  can  then  be  viewed  and  analyzed 
using  other  routines.  The  listing  of  DATRUN  follows  below. 


63 


FILE  DATkuN.bAS 


19  14V 


05 -Uc  t-8  I 


PAGE 


1 


10! 

20! 

30! 

40! 

50! 

60! 

70! 

80! 

90! 


DATKUN.BAS  ....  HAST  DATA  STOKAGE  ROUTINE 

PAUL  BJOKK ,  PROGRAMMER  3-30-81. 

VARIABLE  CHANNEL  INPUT  NUMbEK. 

VAKIABLE  SAMPLE  AND  BUKST  INTERVAL  SET. 
CLOSES  UPON  DISPLAY  TOUCH  OH  AkhAY  OVERFLOW. 
USES  KBYTE  WBYTE  STATEMENT 
IN  3  BYTE  BINARY  INPUT  FORMAT. 

MAXIMUM  SPEED— 200  CHANNELS  PEK  SECOND. 


100! 

110! 

120! 

130! 

140! 

150! 

160! 

170! 

180! 

190! 

200 

210 

220 

230 

240 

250 

260 

2/0 

280! 

290! 

300! 

310 

320 

330 

340 

350 

360 

370 

380 

390 


**  USES  VIRTUAL  akkay  file  ** 

BS ( 0 )  CONTAINS  FILE  FORMAT  INFO.  NO.  OF  CHANNELS. 
bS(l)  CONTAINS  TOTAL  NUMBER  OF  BURSTS. 

A< BURST, CHANNEL)  CONTAINS  VOLTAGE  DATA. 

TB(BUKST)  IS  THE  ELAPSED  RUNNING  TIME. 

**  DVM  HIGH  SPEED  READING  MODE  ** 

REVISION  2...  DVM  IS  SET  TO  THE  2.5  VOLT  RANGE. 

ALTER  LINES  250  AND  870.  P.  BJORK,  8/4/81. 

MULTIPLE  UUhST  SEQUENCES  POSSIBLE  ON  SAME  FILE  8-14-81 

COM  Y 

ESS*CHRS(27)+" I " 

IF  Y= I  GOTO  1740 
□  IM  AD%(9%),  A!  (  10) 

K56-KEY  \  ON  KEY  GOTO  1090 

hF-I.!/8.  !  rF  FOR  20  V  SCALE  IS  !.  RF  FOR  2.5  VOLT  SCALE  IS  1/8 

K4-KEY 

ES$*CHR$(27%)+" I " 


!G0  DIRECTLY  TO  OUTPUT  THE  FILE! 


SET  UP  rt BYTE  COMMAND  STRING  ARhAY. 


AD%<0%)»607% 
AD%(  1 56)  =5  44%+ 1 A 
AD%(2%)=63% 

AD%( 3%) *575% 

AD%<4%)=544%+3% 

AD%(5%)-43% 

AD%(6%)-40% 

AD%(7*)=575% 

AD%( 8%)=576%+ 1  % 


!  ATN,  UNT 
!  ATN,  MTA  DVM 
!  ASCII  TRIGGERS  DVM 
!  ATN,  UNL  DVM 
!  MLA  SCANNER 

!  ASCII  "+"  INCREMENTS  SCANNER 
!  AKBITKAkY  TIME  DELAYER 
!  ATN,  UNL  SCANNER 
!  ATN,  MTA  DVM 


400! 

410  PRINT 
420  PRINT 
430  PRINT 
440  PRINT 
450  PRINT 
460  IF  FI 
4/0  RUN 
480  PRINT 
490  PRINT 
500  PRINT 
510  PRINT 
520  IF 


E5$+"2J" ICPOS<  3  I) 

"DATRUN.BAS  ...’.3  BYTE  2.5  V.  RANGE  HIGH  SPEED  DATA  STORAGE  ROUTINE.1' 


"PLEASE  ENTER  THE  FOLLOWING*" 

"  DATA  FILE  NAME.  ( <CR>  to  view  a  file)  "I  \ 

$<>»"  GOTO  480  ELSE  Y=l  !  TO  VIEW  FILE 


INPUT  FIS 


"  FILE  HEADER.  "I  \  INPUT  FHS 

"  STARTING  CHANNEL,  NUMBER  OF  CHANNELS.  "IN  INPUT  SC%.NC% 

"  BURST  INTERVAL(ms).  "I  \  INPUT  TS 

"  TOTAL  BURST  NUMBER.  "I  \  INPUT  SZS 

SZS»"  THEN  SZ%*300%  ELSE  SZ%=VAL(SZS)  !  MAXImUM  ARRAY  SIZE 


530! 

540!  INITIALIZE  THE  VIRTUAL  DATA  FILE 
550! 
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FILE  DATAUN.BAS 


I5»49  05-Oct-dl 


PACE  2 


5oO  OPEN  FIS  AS  NEW  DIm  PILE  I  SIZE  ( (NCX*  I X)*SZX)/<>4X+2X 
5/0  DIM  #1,  BS(oX)  -  04%,  A(SZ%,NC%) 

550  DIM  AA (NC%) , I A% ( SZ%*3%*NC% ) , TB( SZ%) 

590  NC“NC% 

000  BS  (0X)*NUmS  ( NC )  !  NUMBEk  OP  CHANNELS. 

010  SS(2%)*FI S  !  PILE  NAME 

020  3S < 3%)»FHS  *  PILE  HEADEa. 

o30  dS(4%)*DATES+/  '♦TImES 
040! 

050  PAINT  ESS+'2J' ! ESS+' I  p' ICP0S(2, I > 

000  PAINT  "PAESS  9CAEEN  TO  INITIATE  AEADINGS" 

070  ON  KEY  GOTO  6B0  \  WAIT 

630  K*=KEY  \  AESUME  090 

090  PKINT  ESS+"2J" ICP0S(2, I ) I "SCANNING  CHANNELS  “ISC*I"  TO  " lSC%+NC%-l I M . " 
700  PAINT  CP0S(5,  I  )  I  "DATA  PILE  NAME  IS  "!FI$«".“ 

710  PAINT  CPOS (7, I ) I  "PAESS  SCAEEN  TO  HALT." 

720! 

730!  INITIALIZE  METEA  AND  SCANNEA 

740!  METEA— FILTEA  BYPASS,  A EMOTE  TaIGGEa,  VAAIABLE  SAMPLE 

750!  SCANNEA— SET  LOWEA  AND  UPPEA  I NCAEMENTAL  BOUNDAAIES 

7oO! 

770  IN  IT  POAT  0 
7dO  CLEAa  el 
790  AEMOTE  a  I  e)3 
dOO! 

dIO  IM%=NC%*3%-U  \  UP%=SCX+NC*-I 
820  IN  IT  POAT  0 

d30  PAINT  B3%."*, <$.,)"»  !  AESET  AND  CLEAa  SCANNEA 

840  PAINT  tt3%, SC%I " , ,0"  !  SET  CHANNEL  BOUNDAkIES 

850  PAINT  ti3%,UP%J".,  1“« 

dOO  PAINT  «3%,SCXl","  !  SET  TO  STAATING  CHANNEL 

870  PAINT  ft  1 , "VA2S2F 1T0BQ0! "  !  Al  IS  2.5  V  AANGE.  A2  IS  20  V  AANGE. 

880  MIA  IT  3000 

890  PAINT  ESS+'5m'!  'ACQUISITION  IN  PAOGAESS.'I 
900  INPUT  LINE  »l,01$ 

910! 

920!  PAINT  "ELAPSED  TIME*  0  MINUTES"! 

930!  PAINT  ESS+" 16D" I 
940!  T=TI«E 

950  FOA  J%=0  TO  SZX-2%  !  DO  LOOP  SET  TO  MAX  SIZE 

960  K  3%=3%*NC%* J% 

9/0  ABYTE  (MBYTE  AU%( 0% . .8% ) )  I A%( K3% . .K3X*3%*NC%-I %) »3% 

9dO! 

990  TB(J%)=TIME 

1000!  PAINT  USING  "###*,##", <TB< > )/60 IESS+'7D' I 

1010  ON  KEY  GOTO  1090 

1020  rtAIT  TS-5 

1030  OFF  KEY 

1040  NEXT  J» 

1050  J%=J%-U  !  AEDUCE  J%  POA  NOhMAL  1)0  L(X)P  EXIT 

106o! 

1070!  CONVENT  DATA  TO  aEAL,  CLOSE  FILE,  AND  INFOAM  ON  PILE  STATUS, 
lodo! 

1090  K%=KEY 

1100  IF  K%=o%  GOTO  1200 
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FILE  DATkUN.BAS 
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IIIO  PkINT  ES$*'m'lESS+'2J' 

1120  PkINT  CPOS(3,3)l'kESUM£  DATA  ACQUISITION  ON  TnIS  PILE' 

1130  PkINT  CP0S(5, I ) I'TOUCH  SCkEEN  FOH  CHOICE' 

1140  PkINT  CPOS (7,7)1 'CONVERT  EXISTING  DATA' 

1150  ON  KEY  GOTO  H60  \  WAIT  FOH  KEY 

1160  KX=KEY 

1170  OFF  KEY 

1180  IF  K%>30  GOTO  1190 

1182  PKINT  ESS+'2J'*CPOS<5, I ) I'kESUMING  DATA  ACQUISITION  (5  SEC)' 

1183  WAIT  5000 

1184  PkINT  ES$+'5m'*CPOS<8,5) I ' ACQUI SITION  IN  PKOGhESS' 


1185 

RESUME  1040 

1190 

RESUME  1200 

1200 

JT=J* 

!  JT  IS 

TOTAL  NO.  OF  BUkSTS  MINUS 

1205 

ON  CTHL/C  GOTO  1340 

1210 

T0*TB(UX) 

1220 

PKINT  ESS+'2J'»ES$+'m'»CPoSC3, 

1 >«"CON VERTING 

DATA." 

1225 

CLEAk  «l«3 

1230 

FOH  KX*0%  TO  JT 

1240 

K3X-3%*NCX*KX 

1250 

A(  K%,  NCX)»(T8(KXi)-T0)/l  000 

!  8UKST 

TIME  IS  ON  ThE  NCX  COLUMN. 

1260 

GOSUB  2060 

!  CONVERSION  SUBROUTINE 

1270 

FOk  JX*0  TO  NCX-1 X 

1280 

A(KX.JX)»kA( JX) 

1290 

NEXT  J% 

1300 

PRINT  CHkS (7) 1 

1310  NEXT  K% 

1320  B$<  Ul-NUMSUT+I  )  !  NUMBER  OF  BUkSTS. 

1330  BS(5X)«NUMSUUT,NCX>, “***#♦#. ##">«•"  SECONDS  TOTAL  TIME" 

1340  CLEAk  tf  1  «3 
1350  CLOSE  1 
1360  GOTO  1380 
13/0! 

1380!  SUbkOUTINE  ASK  FOk  INSTRUCTIONS 
1390! 

1400  PkINT  ESS+'2J' lESS+'m' |CP0S<2, I ) I 

1410  PkINT  "KEY  CLOSUkE*  DATA  ACQUISITION  STOPPED." 

1420  PRINT  FIS  *"  IS  FILLED  AFTER  "» 

1430  PkINT  JT+ll"  BURSTS. " 

1440! 

1450!  CHOICE  ENTRY  POINT 
1460! 

14/0  PkINT  CP0S(5,I)«  "CHOOSE  A  FUNCTION." 

1480!  PkINT  ESS+'lm'l 

1490  PkINT  ES$+'5m' ICPOSC o, 6 ) I'EX IT' ICPOSC o, 1 8 ) »' V I  EH' »CPOS( 6,30) » ' AGAI N' 
1500  KX=KEY  \  ON  KEY  GOTO  1510  \  ON  CTRL/C  GOTO  1620  \  WAIT  FOk  KEY 
1510  KX=*KEY  \  RESUME  1520 
la 20  PkINT  ESS+'nr'l 

1530  IF  KX-41  GOTO  lo20  !  THEN  EXIT. 

1540  IF  KX<>49  GOTO  15/0 

1550  PkINT  ESS+'p'j  !  RECYCLE  TO  A  NEW  DATA  FILE. 

1560  kUN 

157u  IF  K4=45  THEN  PkINT  ESS+'p'  ELSE  GOTO  1470 
1580  Y»l  !  VIEW  A  FILE 

1590  KUN 
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1600  PkINT  ESS+'lp" 
lo  10  GOTO  1470 
16  20! 

Io30  PKINT  ESS-fp" 

1640! 

1650  END 
1660! 

1670!  VIKTUAL  PILE  PKINTOUT  SUBkOUTINE 
1680!  A (  I , J )  IS  THE  NUMERIC  AKKAY 

1690!  BS (K  >  IS  THE  PILE  HEADEK  AKKAY 

1700!  -  BS ( 0)  CNTAINS  THE  NUMbEk  OF  COLUMNS  IN  STkING  POkMAT 

1/10!  -  BS<1)  CONTAINS  THE  NUMBEK  OF  HOWS 

1/20!  TB(J)  IS  THE  KUNNING  BUkST  TIME. 

1/30! 

1740  Y=0  \  PKINT  ESS+"p" «ESS+"2J" ICPOS ( 3, 1 ) 

1/50  PKINT  “ENTEK  VIkTUAL  FILE  NAME  FOK  LISTING....  "«  \  INPUT  FJS  \  GOTO  1760 
1760  ON  CThL/C  GOTO  2030  \  ON  KEY  GOTO  2020 
1770  OPEN  FJS  AS  OLU  DIM  FILE  1 

1780  DIM  #1.  B2SI6X)  =  64X,  A5( VAL( B2S ( 1 ) ) , VAL( b2$ ( 0 ) ) ) 

1/90  NC%«VAL<U2S(0)>  \  N6X=VAL( B2S( I ) ) 

I  BOO  PKINT 

1610  PKINT  "FILENAME  "IB2S<2) 

1820  PkINT  B2S ( 3) 

ld30  PKINT  B2  S  (  1 )  l "  BURSTS." 

1840  POK  I%=<4*  TO  5X 
1850  PKINT  82 S ( I % ) 

I860  NEXT  IX 
1870  PKINT 

1880  FOK  I%=0  TO  VAL( B2S (0 ) )- 1 
1890  PKINT  '  CHAN." 1 1 XI 
1900  NEXT  IX 

1910  PKINT  "  TIME  (see.)"  \  PkINT 
1920  FOK  K=0%  TO  NB%-I% 

1930  FOK  J%=0%  TO  NC%- 1 X 

1940  PkINT  NUMS(A5( 1%, J%) , "S###. ####")» 

1950  NEXT  JX 

1960  PKINT  NUMS<A5<I%,NCX>,  "**#*.»") 

1970  NEXT  IX 
1980  PKINT 

1990  PkINT  'PRESS  SCKEEN  TO  CONTINUE." 

2000  ON  KEY  GOTO  2020  \  WAIT 
2010  GOTO  2030 
2020  K*«KEY  \  RESUmE  2030 
2030  CLOSE  1 

2040  PKINT  ESS+"lp"«CP0S(2, I )l"DArKUN"»  \  GOTO  1470 
2o 50  ! 

2060  !  SUbkOUTINE  3BYTE.BAS 

2070  ! 

20d0  !  VEKSION  1.1  12/12/79  J. CHURCH  ILL 

2090  ! 

2100  !  BASIC  VI. 0,  POOS  V1.0 

2110  ! 

2120  !  3 BYTE  WILL  CONVEkT  AN  INTEGER  AKKAY  <  IAX  >  kEAD  IN  PKOM 

2130  !  THE  8502A  IN  HIGH  SPEED  HEADING  MODE  (  3-BYTE  blNAKY  TRANSPEk  ). 

2140  ! 


FILE  UATKUN . BAS 


IV 14V  G5-0ct-d I 


PAGE  5 


2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

22/0 

2280 

2290 

2300 

2310 

2320 

2330 

2340 


CIX-255* 

C2X-256X 

ci«iex/(  iox*kF) 

C2«4096X/< 1QX*HF) 

C3»l04857o/< I0X*hF) 
i 

FOh  IX-K3X  TO  K3X+IMX  STEP  33; 


!  CONSTANT  For  I'S  COMPLEMENT 
!  CONSTANT  FOh  2'S  COMPLEMENT 
!  CONSTANT  FOR  BYTE  I 
!  CONSTANT  FOk  BYTE  2  <2~I2) 

!  CONSTANT  FOk  BYTE  3  <2“20) 


IA%<IX«-2X>*IAX(IX+2X)  ANU  CIS  !  rtEMOVE  EOI  BIT 

IF  IA*(I*>  <  1 2d*  THEN  ISX-U  \  GOTO  2290  !  CHECK  FOk  NEGATIVE  READING 

ISX*-1X  !  SET  SIGN  TO  NEGATIVE 

IAX<IX)«CIX-IAX<U>  !  COMPLEMENT  BYTE  I 

IAXIIX+IXI-CU-IAXC  IX+1XJ  !  COMPLEMENT  BYTE  2 

IAX< IX-*-2X)-C2X-lAX< l**2X)  !  COMPLEMENT  BYTE  3 

IF  lAX(IX)  ANU  32 X  THEN  HA < < I X-<3X)/3X)»OX\GO TO  2320  !  EkkOrt  HAS  OCCUHED 
HA ( ( IX-K3X)/3X)*< IAX< IX )/C I + I AX( IX+ IX ) /C2+ IAX( IX+2X)/C3)*ISX 

i 

NEXT  IX  I  LOOP 

RETURN 

END 
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DPDEV 


In  order  to  examine  the  gyro  performance  over  long-term  drift  runs,  a 
special  acquisition  program  was  developed.  The  routine  acquires  a  burst  of 
multiple  channel  data  similarly  to  DATRUN.  This  burst  is  done  at  high  data 
rate  for  a  short  time  period,  usually  less  than  10  seconds.  These  data  are 
converted  and  statistically  analyzed  for  the  mean  and  sigma  of  each 
channel.  These  values  are  stored  on  file,  and  the  program  returns  to 
acquire  another  burst  of  data.  The  collection  of  means  will  show  the 
low-frequency  drift  character  of  the  gyro,  and  the  sigmas  will  show  the 
high-frequency  character  of  the  gyro  noise.  Therefore,  this  acquisition  and 
real  time  analysis  yields  a  full  bandwidth  examination  of  the  gyro  output 
without  prohibitively  large  data  storage.  The  listing  of  DPDEV  follows. 


PILE  UPUEV.BAS 


19*52  05 -Oct -3 I 


PAGE  I 


10!  DPDEV  ....  EAST  DATA  STOWAGE  KOUTINE 

20!  PAUL  BJOHK,  PKOGhAMMEh  3-30-81. 

30!  VEHSION  3  HE VI SI  ON  i...  DVM  IS  SET  TO  THE  2.5  VOLT  HANGE. 

40!  SThEAMLINED  FOH  'HEAL  TIME'  STATISTICAL  ANALYSIS. 

45!  OUAL  POLAHIZATION  VEHSION  9/2/81 

4/!  DUAL  POLAK IZATION  HEVISION  I  9/10/81  J  HANSE 

50!  P.  BJOHK,  8/4/81. 

GO! 

70  ESS=CHkS(27)+"1“ 

80  J%=4% 

90  DIM  AD%(9%) ,  A1  <  10)  ,MN<J%),TNC  U),SE'(  J*)  .NS(J*),S<  U)  ,ES(  Ul.SI  (  U),S2(  IS) 

100  ON  CTKL/C  GOTO  1240 

110  ON  EHHOh  GOTO  1240 

120  HF  =  I  .  !  HF  FOH  20  V  SCALE  IS  1. 

130  GOSUB  1320  l  3BYTE  INITIALIZATION 
140  STSs'S##. #####' 

150  L 1 = 1 E-5 
loO  KE%*KEY 
170! 

130!  SET  UP  rtdYTE  COMMAND  SThING  AHHAY. 

190! 

200  AD%(0%)=607%  !  ATM,  UNT 

210  AD%<  U)=544%+U  !  ATN,  MTA  DVM 

220  AD%(2%)=63%  !  ASCII  " i "  THIGGEHS  DVM 

230  AD%( 3% 1=575%  !  ATN,  UNL  DVM 

240  AD%(4%)=544%+3%  !  MLA  SCANNErf 

250  AD%( 5% 1=43%  !  ASCII  INCHEMENTS  SCANNEH 

2oO  AD%( 6%1 =40%  !  AhBITHAKY  TIME  DELAYEH 

270  AD%( 7% 1 =575%  !  ATN.  UNL  SCANNEH 

280  AD%(«%1=576%>U  !  ATN,  mTA  DVM 

290! 

300  PKINT  ES$+M2J“ICP0S<3, I H 

310  PHINT  "DPDEV  ....  DATA  ACQUISITION  AND  HEAL  TIME  STATISTICAL  ANALYSIS" 

315  PHINT  '  FOH  OUAL  POLAkIZATION  GYkO  DHIFT  TESTING' 

320  PHINT  CPOS(5,61 I'ENTEH  FILE  NAME  '«  \  INPUT  FIS 

340  PHINT  CPOS( 7,1211'FILE  HEAUEH  'I  \  INPUT  FHS 

350  PHINT  CPUS! 8,121 !' STAKTI NG  CHANNEL  '»  \  INPUT  SC% 

360  PHINT  CP0SI9, 12) I'The  input  channels  are  assumed  to  be  '»SC%;'  and  'ISC%+1% 
370  PHINT  CPOS<  10,  12U'BUHSTS  PEH  CYCLE  'l\  INPUT  NB% 

380  PHINT  CPOSf  11,12) I 'NUmBEH  OF  CYCLES  'l\  INPUT  CY* 

390! 

500!  INITIALIZE  THE  VIHTUAL  DATA  FILE 
510! 

515  NC%=2% 

520  OPEN  FIS  AS  NEW  DIM  FILE  I  SIZE  ( 5%*CY%)/64%+2% 

530  DIM  #1,  BS(6%)  =  64%,  A<CY%,4%) 

540  DIM  K A< NC%) . I A%< NB%*3%*NC%1 

550  BS ( 0% 1 =NUMS ( 4%)  !  NUMbEH  OF  CHANNELS  STOhED  ON  FILE. 

5oO  BS(2%)=FIS 

570  BS ( 3% )=FHS 

530  BS(4%)=UATES+"  m+TIMES 

590! 

oOU!  INITIALIZE  METEH  AND  SCANNEH 

610!  mETEH— FILTEH  BYPASS,  HEMOTE  TKIGGEH,  VAHIABLE  SAwPLE 

o20!  SCANNtH— SET  LOAEK  AND  UPPEH  INCKEMENTAL  BOUNDAH I ES 
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FILE  DPDEV.8AS 


19*52  05-Oct-d I 


FACE  2 


030! 

040  IN  IT  worn  0 
050  CLEAh  It  I 
000  hEMOTE  tl«3 
0/0  IMX»NCX*3X-IX 

030  UPX=SCX+NCX-I  !  UPPEk  bOUND 

090  INIT  POkT  0 

700  PHINT  m3X,"*,<S,,)»i  ! 

710  PHINT  83X,SCXI',,,0“  ! 

720  PHINT  83*. UP**",, l“l 

730  PHINT  a3X,SCXl",  "  ! 

740  PHINT  tt I  ,»VH2S0F1T0BQ0!M  ! 

750  INPUT  LINE  81,01 S 

700  PHINT  CP0SC13,  HI"  TIME  (MIN) 

7o5  PHINT  CPOSC 14, 19) l"MEAN  SIGMA 

770  PHINT 
780  TO=TIm£ 

790! 

800!  CYCLE  ACOUISITION  AN!)  ANALYSIS  L(X)P 
310! 

d20  ON  KEY  GOTO  1170 
330  fOk  LX-OX  TO  CYX-IX 
840  ACL%.OX)*<TIME-TO)/< 1000*00)  ! 

845  FOk  JX=OX  TO  IX 
850  S(JX)=0. 

300  SS(JX)=0. 

805  NEXT  JX 

870  OFF  KEY 

880  FOH  KX»OX  TO  NBX-IX 
890  K3X=3%*NCX*KX 

900  HBYTE  (MBYTE  ADX(OX..oX)>  I AX(K3X. .K3X+3X*NCX-1 X) *3X 
910  NEXT  KX 
920  ON  KEY  GOTO  1170 
930  FOH  KX=OX  TO  NBX-IX 

940 

950  GOSUB  1410  !  CONVEKSION  SUBHOUTINE 

900!  FOH  JX=  1  TO  NCX 

970!  NS( JX )  =  (HA( JX-I X) -MN<  JX ) )/SF ( JX)  !  SINUSOIDS 

930!  NEXT  JX 

1015  FOH  JX=OX  TO  IX 

1020  S(JX)=S< JX)+HA< JX) 

I U40  SS( JX )  =  SS( JX)  +  kA( JX)*kA<  JX) 

1045  NEXT  JX 

1050  NEXT  KX 

1053  PkINT  USING  "###.##", A< LX , OX) , 

1055  FOH  JX*OX  TO  IX 

1080  SI(JX)=S(JX)/NbX 

1070  S2( JX ) =SQH ( ( SS( JX )/NBX ) -SI ( JX)*2) 

1080  A( LX, 1 X*2*JX) «SI ( JX )  !  1,3  AkE  MEAN  OF  DET  I  AND  2 

1090  A(LX,2X+2*JX)-S2(JX)  !  2,4  AKE  SIGMA  OF  UET  I  AND  2 

1100  PhINT  USINO  STS.SI  <J-»)I"  "IS2(JX), 

1103  NEXT  JX 
1105  PhINT 

1110  IF  KEXoO*  GOTO  1220 
1120  NEXT  U 


HESET  AND  CLEAH  SCANNEH 
SET  CHANNEL  BOUNDAHlES 

SET  TO  STAHTING  CHANNEL 

Hi  IS  2.5  V  kANGE.  H2  IS  20  V  HANGE. 

DET.  1 ' l TAB (65X1  I "DET.  2" 

MEAN  SIGMA" 


TIME  STOHED  IN  MINUTES 


FILE  DPUEV.BAS 


19 «52  05-Oct-BI 


PACE  3 


1130! 
1140 
1150 
1 160! 
II  70 
1180 
1190 
1200! 
1210! 
1220 
1230 
1240 
1250 
1260 
1270 
1  280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1260 
13  70 
1300! 
1390! 
1400! 
1410 
1420! 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
15  20 
1530 
1540 


L*«LX-U 
GOTO  1210 


!  FOH  NORMAL  DO  L(X)P  EXIT 


PRINT  "KEY  CLOSURE <  ANALYSIS  .*ILL  CONTINUE  TO  FINISH  CYCLE.' 

KEX-KEY 

RESUME 

CLOSE  FILES  AND  END 

B$< IX)*NUMS<LX)  !  NUMBER  OF  CYCLES 

BS<5%)*NUMS(A<LX,0X))+'  MINUTES  TOTAL  TIKE' 

CLEAN  til  1 03 
CLOSE  I 

IF  EKko  0  THEN  PKINT  'EKHOK  ' lEKKI'  OCCUKEU  ON  LINE  'JEKL  \  GOTO  1270 

CLOSE  1 

END 

i 


!  SUUHOUTINE  3BYTE  INITIALIZATION. 

{ 

CIX-255X 
C2*»256X 
C I-I6X/I IOX*KF) 

C2=»4096X/<  10X*KF) 

C3= 1 048576/1 10X*KF) 

KETUKN 

SUBROUTINE  3EYTE  CONVERSION 
FOK  I X=K3X  TO  K3X+IM*  STEP  3% 


!  CONSTANT  FOK  I 'S  COMPLEMENT 
!  CONSTANT  FOK  2'S  COMPLEMENT 
!  CONSTANT  FOK  BYTE  1 
!  CONSTANT  FOK  BYTE  2  (2“I2> 

!  CONSTANT  FOK  BYTE  3  (2“2 0) 


IAX<I*+2%}=IA%<I*+2X)  AND  Cl* 

IF  I A%( IX )  <  128X  THEN  ISX=IX  \  GOTO  1490 
IS*=-I* 

IAX( I*)=C 1X-IA*( IX) 

IAX(IX+I*)=C1*-IAS(  IX+IX) 

IAX(IX+2*)=C2*-IAX( IX+2X) 

IF  I  A* ( I  * )  AND  32*  THEN  K  A ( ( I X-K3  X )/3X) =OX\GOTO  1520  !  EHKOK  HAS  OCCUKED 
KA ( ( I X-K 3* ) /3X )  =  ( I  A* ( IX)/C !♦ IAX ( IX+ 1X)/C2+ IA*( IX+2*)/C3)*ISX 


KEMOVE  EOI  BIT 
CHECK  FOK  NEGATIVE  READING 
SET  SIGN  TO  NEGATIVE 
COMPLEMENT  BYTE  I 
COMPLEMENT  BYTE  2 
COMPLEMENT  BYTE  3 


NEXT  IX 

RETURN 

END 


!  LCX1P 
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DPSTAT 


This  is  a  simple  statistical  analysis  package.  The  data  is  input  from  the 
diskette.  The  routine  can  accept  data  from  either  DATRUN  or  DPDEV.  The 
routine  spools  through  the  data  and  calculates  the  mean,  sigma,  maximum 
value  and  minimum  value  of  each  channel  selected  by  the  operator.  These 
results  are  output  to  the  printer.  The  routine  then  goes  back  and 
calculates  a  best  fit  straight  line  for  the  data  using  time  as  the 
independent  variable,  removes  the  line  from  the  data,  and  calculates  the 
residue.  The  residue,  slope,  and  intercept  are  output  to  the  printer.  The 
units  of  the  channels  are  preserved,  so  if  the  data  were  stored  as  volts 
every  10  seconds,  the  slope  would  be  in  volts  per  second  and  all  the  other 
outputs  would  be  in  volts 


FILE  DPSTAT.  BAS 


I V  *3  I  05 -Oct-8 I 


PAGE 


I 


10!  OPSTAT . . ..DUAL  POLARIZATION  STATISTICAL  DATA  ANALYSIS  ROUTINE. 

20!  PAUL  dJOHK,  PROGRAMMER. 

30!  VERSION  I,  REVISION  3,  8/5/81. 

40!  DUAL  POLARIZATION  VERSION  9/10/81  JOEL  HANSE 

50!  FUNCTIONS! 

60!  I.  CALCULATES  MEAN  AND  SIGMA  ON  CHANNEL  IN  A  OATRUN  TYPE  FILE 

70!  2.  REMOVES  A  LEAST  SQUARES  STRAIGHT  LINE. 

80!  3.  CALCULATES  MEAN  AND  SIGMA  OF  RESIDUAL. 

90! 

100  E3S=CHKS ( 27 )*“ I " 

105  CHS=ESS+'ri'+ESS+'2J'  !  CLEAR  SCREEN  AND  HOME 

106  LFS«CHKS(|0)  !  LINE  FEED 

107  FFS*CrihS( 12)  !  FORM  FEED 

1 08  OPEN  'KBOl'  AS  NEW  FILE  4 

109  OPEN  'KBH'  AS  NEW  FILE  5 

110  PRINT  CHS ICP0SI3, 1 ) » 

120  on  Error  goto  1030 
130  ON  CTHL/C  goto  1030 
140  J%*5% 

150  DIM  S  (  JX  ) ,  SSI  JX ) ,NC%( JX), YXIJX) . YM I  JX) , SX( JX) . XX ( JX) . XY ( JX) ,M ( JX ) . B( JX) 

loO  STS*' S*#. #####' 

170! 

180  PRINT  ESS+'2J'lCPOS(3, 1)1 

190  PRINT  'DPSTAT ... .DUAL  POLARIZATION  STATISTICAL  DATA  ANALYSIS  ROUTINE.' 
200  PRINT 

210  PRINT  '  This  routine  performs  three  functions' 

220  PRINT  '  I.  Calculates  mean  and  sigma  of  any  channel  of  a' 

230  PRINT  '  Standard  (DATHUN)  type  virtual  data  file' 

240  PRINT  '  2.  Removes  least  squares  straight  line' 

250  PRINT  '  3.  Calculates  mean  and  sigma  of  residual' 

260  PRINT  '  Items  2  and  3  can  be  disabled  on  command' 

270  PRINT  CPOSI I4,5)«'ENTER  DATA  FILE  NAME.  'I  \  INPUT  FIS 
280  OPEN  FIS  AS  OLD  DIM  FILE  1 

290  DIM  #I,BSI6X)*64X,A(VAL(BS( IX)), VALIBSIOX))) 

300  NBX=VAL(BS( IX)) 

320  PRINT  CHSICP0S(5, I ) I'FILt  ' * H I S * '  CONTAINS  'INBX-IX!'  DATA  POINTS' 

330  PRINT  CPOSI 7 ,1)1 'FOR  PARTIAL  FILE  DATA  ANALYSIS  ENTER' 

340  PR INTCPOSI 8. 5 > I 'LINE  NUMBER  OF  FIRST  DATA  POINT  TO  CONSIDER  ' 

350  PR INTCPOSI 9 , 5 ) I 'AND  NUMBER  OF  POINTS  TO  ANALYZE  '» 

360  INPUT  FPX,N 

3o5  IF  FPX+N>NbX  GOTO  330 

367  LPX=FPX+N-I 

370  Print  CPOSI II, I ) I' ENTER  X  AXIS  CHANNEL  For  LEAST  SQUARES  LINE  REMOVAL  ' l 
375  PRINT  CPOSI 12,5) »'C/R  IMPLIES^NO  LEAST  SQUARES  LINE  REMOVAL  'I 

380  INPUT  XCS 
385  PRINT  CHS! 

390  IF  XCS<>"  THEN  XCX=VAL(XCS>  \  GOTO  420 
400  XCX*-IX 

410  PrtINT  CPOSI 7, 1 ) I 'CALCULATING  ONLY  RAW  DATA  STATISTICS.' 

420  PRINT  CP0SI8.5) »'  ENTER  COLUMNS  TO  BE  ANALYZED  '»  \  INPUT  NCS 
430  CTX=LEN (NCS ) 

440  FOR  IX*  IX  TO  CTX 

450  NCX(U)=VAL(MID(NC$,IX,  U>) 

460  NEXT  IX 

4/0  GOSUB  1400  !  ZERO  ARRAYS 
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HI L£  DPSTAT.BAS  I9»3l  05-Oct-8I  PAGE  2 

49U  PrtINT 

500  PrtINT  CP05(  I  1,5)  « ‘'CALCULATING  MEAN  AND  SIGMA  ON  FILE  '|FIS 
510  FOH  IX=FPX  TO  LPX 
520  FOrt  JX=IX  TO  CTX 
530  CNX-NCXCJX) 

540  Y=A<IX,CNX) 

550  GOSUB  1250  !  CALCULATE  SUMS 

500  NEXT  JX 
570  NEXT  1% 

580! 

582  PrtINT  CHS I 

584  HOk  10*  4  TO  5 

530  PhINTmIO.LFSJLFSILFSJLFSJLFS 

590  PHINTMIO,'  STATISTICS  ON  DATA  FILE  '«FIS|'.' »LFSlLFSlLFS 

OI0  PrtINTMIO, 'COLUMN  NO.  MEAN  SIGMA  MINIMUM  MAXIMUM' 

020  PH  I NT# 1 0 

630  FLX=IX  !  FLAG*  PrtINT  MAX  AND  MIN. 

640  GOSUB  1100  !  PH  I  NT  rtESULTS. 

642  NEXT  10 

650  IF  XCX*-IX  GOTO  1010  !  STOP  IF  NO  STKAIGHT  LINE  FIT  DESIHED. 

660! 

070!  CALCULATE  STHAIGHT  LINE  LEAST  SQUArtES  FIT. 
o80! 

685  GOSUB  1400  !  ZEKO  ARRAYS  AGAIN  FOR  STHAIGHT  LINE  FIT 

690  PRINT  CHSICPOSI 11,5)1 'CALCULATING  STHAIGHT  LINE  FIT.' 

700  FOH  IX=FP%  TO  LPX 


710 

For  JX*IX  TO  CTX 

720 

CNX=NCX( JX) 

730 

X-A(IX.XCX) 

!  X  AXIS  CHANNEL. 

740 

Y=A( IX.CNX) 

750 

GOSUB  1250 

!  CALCULATE  SUMS. 

760 

GOSUb  1330 

!  CALCULATE  LINE  SUMS. 

770 

NEXT  JX 

780  NEXT  IX 

790  Fok  JX»IX  TO  CTX 

dOO  M(JX)=<XY<Jx)-SX<JX)*S<JX)/N)/(XX<JX)-<SX(JX)~2>/N)  !  SLOPE. 

dIO  8(JX)*(S(JX)-M(JX)*SX(JX))/N  !  Y  INTERCEPT. 

820  NEXT  JX 

d30  PrtINT 

840! 

850!  CALCULATE  STATISTICS  ON  RESIDUE. 

800 ! 

870  GOSUB  1400  !  ZEHO  AHHAYS 

ddO  Fok  IX*FPX  To  LPX 
890  FOH  JX*1X  TO  CTX 

900  CNX*NC%( JX) 

910  Y=A<IX,CNX)-(m(JX)*A< IX, XCX )  +  B<  JX ) )  !  REMOVE  LINE  FROM  Y  POINT 

920  GOSUB  1250  !  CALCULATE  SUMS. 

930  NEXT  JX 

940  NEXT  IX 

950  PRINT  CHS  I 

952  Fok  10  *  4  To  5 

954  PRINT  9I0,LFS|LFS»LFS»LFS 

9o0  Pk I  NT# 10,  '  RESIDUE  AND  EXTRACTED  LEAST  SOUAHES  LINE  STATISTICS.' 

970  HS= 'COLUMN  NO.  MEAN  SIGMA  SLOPE  Y  INTERCEPT' 
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97 2  PR  I  NT  #10, HS 
9ao  gosuo  lloo 

990  PR I NT»IO,  'CALCULATIONS  COMPLETED  on  FILE  '»FISI'.' 

992  NEXT  10 
IOOC! 

lulO!  CLOSE  FILES  AND  END  uF  PkOGrtAM 
1020! 

1030  IF  EkROO  THEN  PRINT  'EKROk  '  «E(<r I 'OCCURED  ON  LINE  'iEhL*'.' 

1040  CLOSE  1 ,4,5 
I05u  END 

I  OoO ! 

lo7u! 

1030!  CALCULATE  wEAN  AND  SIGMA,  PRINTOUT  ROUTINE. 

1090! 

I I  00  Fort  J%=IS  TO  CTs 

1110  CN%=NC*(JA) 

1120  SI =S< JX)/N 

1130  S2=SQk<<  SS(  J*)/N)-Sr2) 

1140  PRINT  #10 

1150  PRINT  #IO,CNA»"  "I 

IloO  PrINT  #10, USING  STS, SI  lTAB<  26 ) « S2l  TAlil  40 ) « 

1170  Ir  FU=U  THEN  Prt  I  NT#  10,  USING  STS  ,  YmCCN*  > *  T  Ab  ( 54  ) ;  YX  (CNA ) \GOTO  1190 
lldO  PrtINT»IO,  USING  STS,  4<  JS)  «TAB<  54)  I  bUX) 

1190  NEXT  J* 

1200  PrINT#IO 
1210  RETURN 
1220! 

1230!  SUBROUTINE  CALCULATE  SUMS. 

1240! 

1250  3(J*)*S(J;4H-Y 

1 2oo  ss<ja)=*ssu*)-«-y~2 

1270  IF  YX  (CN4)  <Y  THEN  YX(CN'A)=Y 

1230  IF  YM(CN/®)> Y  THEN  Y'.,<CN*)=Y 

1290  rETUrN 
1300! 

1310!  SUbkOUTINE  CALCULATE  LINE  SUMS. 

1320! 

I23U  SX<J*)=SX<JA)+X  !  X  SUM 

1340  XY(JS)=XY< JA)+X*Y  !  XY  PRODUCT  SUM 

1350  XX(J*)=XX< JX)+X~2  !  X  AXIS  SUM  Or  SOUARES. 

1300  RETURN 
1370! 

1330!  SUbrtOUTI NE  ARRAY  ZEROING  AND  INITIALIZATION 
1390! 

1400  FOR  JX>=U  To  CTA 
1410  FLX=Oa 

1420  3UA)=0.  !  INITIALIZATION 

1430  SSI JA) =0 • 

1440  3XIJA)=0. 

1450  XXIJA)=0. 

1460  XY I J A) =0. 

14 7o  ONa=NCA(JA) 

1430  YX(C:U)=A(UA,CNM) 

1490  YM(CNA)=A(OA,CIU) 

IdOO  TEXT  JA 


!  SUM. 

!  SUm  OF  SQUARES. 

!  MAXIMUM  VALUE. 
!  MINIMUM  VALUE. 


FILE  OPSTAT.oAS 


19*31  05-0ct-0l 


PACE  4 


1 5 1 o  re  Turn 

1520  END 


A. 

0 


i 


76 


RATMAX 


This  simple  routine  spools  through  a  data  file  that  has  been  generated  by 
DATRUN  and  determines  the  maximum  and  minimum  values  for  each  channel.  The 
results  are  printed  out  and  also  stored  on  diskette  for  later  use.  This 
routine  is  particularly  useful  in  determining  the  intensity  scale  factor  of 
the  gyro. 

Program  RATMAX  will  normalize  any  dual  polarization  data  file  stored  as  two 
single  polarization  detector  outputs  to  the  dual  polarization  mode  12-11, 
and  store  the  normalized  data  with  time  on  an  output  file  with  a  .NOR 
extension.  Full  rate  scale  is  normalized  to  the  sinusoidal  -1  to  +1  for 
easy  computational  reference. 
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10!  HATMAX . .  . HOUTINE  PACKAGE  TO  FIND  MAX  AND  MIN  OF  FIbEH  CM)  HATE  ANGLE  FILES 
26!  HEVISION  I,  PkOGkAMMEH  PAUL  BJOkK,  10-2/81. 

30!  SUBROUTINES* 

40!  LIMITS. ...HATE  HUN  SCALING  LIMIT  CALCULATION. 

70!  MULTIDETECTOH  DATA  INPUT  FILE  FOHMAT— 'DATkUN' ,  EXTENSION  '.DAT'. 

90 ! 


100! 

INPUT.DAT 

OUTPUT. ANG 

CHANNEL  DESCHIPTION 

110! 

0 

0 

Kate  or  ether. 

120! 

2 

1 

Herizental  Pelarlzatlen  angle. 

130! 

1 

2 

Vertical  Pelarlzatlen  angle. 

140! 

Last 

3 

Elapsed  time 

150! 

160  ESS=CHhS ( 27)+' I ' 

165  ON  CTHL/C  GOTO  510 

1  70  DIM  F(8X, 1 X) , SF( 4X) , 

mN(4X),mX(4X) 

175  PkINT  ESS+'2J' 

160  PHINT  CP0SI3X, IX)I'HATMAX....KATE  ANGLE  SCALING  LIMITS  HOUTINE.' 

190  PHINT  CP0SI5, 5) I'ENTEH  DATHUN  HATE  DATA  FILE  NAME' I 

200  INPUT  FIS 

210  YX=INSTH< IX, FIS,'.') 

220  F0S=LEFT(FIS,YX-1X)+'.N0k' 

225  PHINT  CP0SI6, 8) I'ENTEH  FIHST  DETECTOH  COLUMN  NUMBEH  (HATE<=1 ,DhIFT«0)  '» 

227  INPUT  SCX 

230  PriINT  CPOSt  7 ,5 ) I'ENTEH  DATHUN  LIMITS  FILE  NAME' I 

240  INPUT  FHS 

245  IF  FHS="  THEN  Fk$=FI$ 

250  YX=INSTH( IX, FHS,'.') 

26 0  IF  YX=OX  THEN  YX=LEN(FHS )+ I  % 

270  FH$=LEFT(FH$,YX-IX) 

280  ON  EkHOK  GOTO  430 

290  OPEN  FHS+'.MXM'  AS  OLD  DIM  FILE  6  !  TEST  FOH  LIMITS  FILE. 

300  DIM  96, FFS(6X)=64X,F0< I0X.5X)  !  JUMP  TO  430  IF  NO  EXISTING  LIMITS  FILE. 
310!  ON  EHHOH  GOTO  510 

320  GOSUB  4000  !  NOkMALIZATION 

360! 

370!  PHOGHAM  COMPLETED*  riOUSECLEAN. 

360! 

390  CLOSE  I, 2. 5, 6 
395  PHINT  CPOS( I 2X  5X)I 

400  PHINT  'KATE  ANGLE  FILE  'iFOSI'  HAS  BEEN  FILLED  AND  CLOSED.' 

410  END 
420! 

430!  LIMITS  FILE  EHHOH  POINT 
440! 

450  IF  EHK  <>  305  GOTO  510  !  EHHOH  305  IS  NO  EXISTING  FILE  ON  DEVICE 

4o0  kESUME  465 

465  PHINT  CPOS(7X,5X)l 

4od  PkINT  'NO  AVAILABLE  LIMITS  F I Lh . ..CALCULATING  NEW  SCALING  kATIOS.' 

470  GOSUb  2000  !  CALCULATE  NEW  LIMITS  i-ILE. 

4dO  GOTO  540 
500! 

510!  EHHOH  NOTIFICATION 
520! 

530  PHINT  'EHHOH  NUMBEk  'iEHHI'  HAS  OCCUkED  ON  LINE  'lEnLl'.' 

540  CLOSE  1 ,2,5,6 
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550  END 
5o0! 

2000!  SUBROUTINE  LIMITS. ...DETERMINES  SCALING  LIMITS  FROM  A 
2010!  HATE  NUN  FOR  USE  B Y  THE  'ANGLE'  SUBROUTINE. 

2020!  INPUT  FILENAME  IS  IN  'OATRUN'  FORMAT  WITH  A  '.DAT'  EXTENSION. 

2030!  OUTPUT  FILE  IS  'Katerun  fl lename.MXM'  ON  THE  SAME  FLOPPY  DISC. 

2130! 

2145  CN-O.  !  .2  !  FILTER  CONSTANT. 

2150  OPEN  !■"«<*♦' .OAT'  AS  OLD  DIM  FILE  5 

2160  DIM  #5,85$<6X)=64X,A5(VAL< B5$< IX ) ) , VALC b5$<0X ) ) ) 

2170  Fh$«FR$+'.MXm' 

2 1 ao  OPEN  Fk$  AS  NEW  DIM  FILE  6 

2190  DIM  #6,FF$(6X)=64X,FD(9Xt5X)  !9  OR  10? 

2200  NbX«VAL(B5$< 1 )) 

2210  NCX*VAL<B5$<0)) 

2240  For  JX=SCX  TO  IX+SCX 

2250  TI=A5(0X,JX)  !  INITIALIZE  MAX  AND  MIN. 

2260  MX(JX)=*T1 
2270  MN<JX)*TI 

2290  NEXT  JX 

2300!  FILTER  AND  LIMITS  CALCULATION  LOOP. 

2310! 

2320  FOR  IX=OX  TO  NBX-IX 
2330  FOR  JX-SCX  TO  IX+SCX 

2340  FIJX, IX)=A5(IX,JX) 

2360  IF  MX<JXXF<JX,  IX)  THEN  MX< JX)=F< JX. I  X)  !  TEST  FOR  MAXIMUM. 

2370  IF  MN(JX)>F(JX, IX)  THEN  MN< JX)=F(J%. IX)  !  TEST  FOH  MINIMUM. 

23  aO  NEXT  JX 

2390  NEXT  IX 
2400! 

2410!  WRITE  DATA  TO  THE  OUTPUT  FILE. 

2420! 

2425  PRINT  CPOS(  IIX, IOX)l  '  CHANNEL  MIN  MAX' 

2430  FOR  JX=SCX  TO  SCX+IX 
2440  FD(JX,  IX)=MN(JX) 

2450  FD( JX, 2X)=MX( JX) 

2455  PRINT  , JX.MNI JX) ,MX< JX)  !  KILL  THIS 

2460  NEXT  JX 

2470  FD ( 9X , OX ) = A5 ( N  BX , NCX ) 

24 aO  FDI9X, I X)=A5< NbX, NCX )/NfaX  !  THIS  MAY  BE  ADJUSTED. 

2490  FD(9X,2X)=NBX 

2520  FF$(0X)*B5$(2X) 

2530  FFS  (  IX)  =85$  ( 3x ) 

2540  PRINT  CPOS(aX, 5X) I'ENTER  TWO  LINES  OF  INFORMATION  DETAILING  LIMITS.' 
2550  PRINT  CP0S(9X,  10X)  I  \  INPUT  F'SUX) 

2560  Print  CPOS< I0X, IOX)«  \  INPUT  FF$C2X) 

2570  PRINT  CP0SI7X, IX) «ES$+'2K' I 

25 BO  PRINT  CP0SC7X,5XM  'LIMITS  FILE  MFtm'  IS  COMPLETE.' 

2585  PRINT  CPOS(8X,5X) «ES$+'J'I 
2590  CLOSE  5 
2600  RETURN 
26 10  END 

4000!  SUBROUTINE  NORMAL ....  RATE  SCALING  LIMITS  DATA  NORMALIZATION. 

4010!  SCALES  AND  NORMALIZES  EACH  DATA  CHANNEL  TO  UNITY  USING  MAX  AND 
4020!  MIN  DATA  STOktD  IN  FR$='Raterun  fi lename.MXM' . 
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4060!  OPEN  DATA  FILES.  LIMITS  FILE  *6  ALPEAOY  OPEN. 

4070! 

4080  OPEN  FIS  AS  OLU  DIM  FILE  I  !  INPUT  FILE. 

4u90  DIM  #1 ,BIS<6)«64*,AI<VAL(BISC I ) ) , VAL( 81  $ (O ) ) ) 

4093  NC%»VAL<BIS<0)) 

4095  NBX*VAL(BI$< I )) 

4100  OPEN  FOS  AS  NEW  DIM  FILE  2  SIZE  2*NdX/64X+2X  !  OUTPUT  FILE. 

4110  DIM  #2, B2$( 6X) *64X, A2( VAL<  B I  $ ( I  ) ), IX) 

4  120! 

4130!  FILTEH  INITIALIZATION 
4140! 

4200  DIM  AF( NBX, 2%) 

4210  FOH  JX=IX  TO  2X 

4220  SF(JX)»<FDUX,2X)-FD(JX,  1X>)  !  SCALE  FACTOP  CALCULATION. 

4240  MN(JX)=FD(JX, IX)  !  MINIMUM  FOP  EACH  CHANNEL. 

4250  NEXT  JX 
4260! 

4270!  FILTEK  AND  NOrtMALIZATION  PPOCESSING  L(X)P. 

4280! 

4290  FOP  IX=OX  TO  NBX- IX 

4330  AF(IX,OX)»<AIUX.SCX)-MN<  U))/SF(  IX)-. 5  !  SCALE  DOWN. 

4332  AF(IX,1X)*<A1<IX,SCX+IX)-MN(2X))/SF(2X)-.5  !  SCALE  DOWN. 

4345  NEXT  IX 

4350  FOH  IX=OX  TO  NBX- IX 

4353  IF  SCX»1X  ThEN  A2(  IX,0X)*AM  IX.OX)  S  GOTO  4380  !  PATE  DATA  ON  LINE  0 
4355  A2<IX,0%)=AI (IX.NCX) 

4380  A2(I%,  IX)»AF(U,  l  X)-AF(  U.OX) 

4395  NEXT  IX 
4401  ! 

4403  62S(0X)=NUMS< U) 

4404  82 $( IX) =81 S ( IX) 

4405  82$<2X)-F0$ 

4406  82  S ( I  X) «BI $ ( IX) 

4407  FOH  JX=3X  TO  6X 

4408  82 $(JX)-B I $(JX) 

4409  NEXT  JX 
4410! 

4450  CLOSE  I ,2,5.6 
4460  PETUHN 
4470! 


HPL0T1 


HPLOTl  is  a  general-purpose  plot  routine  that  drives  an  Hp  7225A  plotter. 

The  input  data  format  is  the  same  as  the  output  of  DATRUN.  The  operator  has 
control  over  the  axis  scaling,  titles,  and  which  data  to  plot.  The  routine 
is  capable  of  plotting  up  to  4  channels  of  data  on  a  single  graph. 
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10  !  'HPLOTI'  ....  ROUTINE  TO  00  MULTIPLE  PLOTS  ON  HP  7225A 

11  !  JOEL  HANSE  6/10/80 

12  ! 

13  !  PLoTTEk  IS  CONNECTED  I EEE-4bd  PONT  '105' 

15  !  DATA  IS  ASSUMED  TO  bE  A  DATHUN  FORMAT  VETUAL  ARRAY 
17  ! 

Id  ESS*CHRS<  27 )♦' I ' 

19  PRINT  ESS+'2J'lCPOS(3, I ) 

20  ON  CTkL/C  GOTO  1240 
30  ON  ERROR  GOTO  1240 
32  I  NIT  PO«T  I 

35  PHINT  'HPLOTI  ....HP  7225A  PLOTTING  ROUTINE  FOR  DAT  RUN  FORMAT  FILES.' 

37  PRINT 

40  PRINT  'ENTER  FILE  NAME' I \  INPUT  FIS 

50  PC%=I  !  PLOT  COUNTER  FOh  MULTIPLE  PLOTS 

60  PRINT  FIS  I'  IS  ASSUMED  TO  bE  A  UATRUN  VIkTUAL  ARRAY  FILE.' 

05  PRINT  'ENTER  PLOT  HEADING' I  \  INPUT  PHS 
7o  DIM  X  <500 ) , Y  <  500 ) 

80  OPEN  FIS  AS  OLD  DIM  FILE  I 

90  DIM  §  1 , b2S (6 )=64%, V <VAL<b2S(l  )),VAL(B2S<0))> 

100  IF  PCX<> I  GOTO  140 

110  PRINT  'ENTER  COLUMN  FOR  X  AXIS'»\INPUT  IXX 
120  PRINT  'IS  X  AXIS  ThE  TACh  OUTPUT' I \INPUT  KRS 
130  IF  LEFT ( KRS , I )  =  'Y'  THEN  KR=60/1 . 10  ELSE  KH=1 
140  PRINT  'ENTER  COLUMN  FOR  Y  AXIS'JMNPUT  IY% 

150  NCX=VAL( B2S (01) 

100  IF  I Y%<0  THEN  1240 

170  IF  IXX>NC*  OR  IY%>NC*  THEN  100 

172  IF  PC*<»1  GOTO  240 

180  NP2=VAL(B2S<  ).)) 

190  PRINT  'DATA  POINTS  IN  FIL£«'INPX 

200  PRINT  'ENTER  *  OF  DATA  POINTS  TO  PLOT'*\INPUT  NPX 

210  !  READ  IN  THE  X  AND  Y  AXIS  AND  COMPUTE  MAX  AND  MIN  VALUES 

230  XX*-IE30\Xm=IE30 

240  YY=*-1  E30\YM=  1 E30 

241  V  ( NP%,  IYX)=V(NPX- 1%,  I  YX 1 
250  FOR  JX=  1  TO  NPX- 1  % 

260  IF  PCX  <>  IX  GOTO  300 

270  X(JX>)=V(J%,IXX)*KH 

280  IF  X(JX)>XX  THEN  XX=X<JX) 

290  IF  XUXXXM  THEN  Xrt=X(JX) 

300  D3=ABS<VCJ%+IX,m>-V<JX-U,IYX>>  !  WILD  POINT  EDITOR 
302  D2=ABS(  V(  JSS+  l%,IYX)-V(JX,IYX)> 

304  D1»ABS(V(JX,  IYX > -V < J%- IX,  IYX  >> 

306  IF  DI>50*D3  AND  D2>5 0*U3  THEN  Y< JX >=V< JX- I %, I Y*D  ELSE  Y<  J*)»=V<  J%,  IYX) 
310  IF  Y ( JX )> YY  THEN  YY”Y <  JX) 

320  IF  Y ( JX) < YM  THEN  YM=Y(Jx) 

330  NEXT  JX 

335  Y ( OX ) =V (UX, I YX) \X (OX ) =V (OX, IXX)*KR 

340  IF  PCX<>1  GOTO  370 

350  PRINT  'X  MAX  AND  MIN  VALUES  ARE' I  XX, XM 

Js>0  PRINT  'ENTER  X  AXIS  mAX  AND  «IN  FOR  THE  PLOT'lMNPUT  SS.SM 

361  IF  SS>Sm  THEN  XX«SS  ELSE  XX=Sm 

Jo2  IF  SS>SM  Then  XM*SM  ELSE  Xm»SS 

370  PRINT  'Y  MAX  AND  ,.iIN  VALUES  ARE'  I  YY,YM 
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3dO  PRINT  'ENTER  Y  AXIS  MAX  AND  MIN  POP  THE  PLOT'IMNPUT  SS.SM 

381  IP  SS»SM  THEN  YY  =  SS  ELSE  YY=SM 

332  IP  SS>SM  THEN  YM-SM  ELSE  YM=SS 

390  IP  PCXolX  GOTO  420 

400  XL=6.\YL=6. 

410  PKINT  'ENTER  X  AXIS  LA EEL' \ INPUT  LXS 
420  PRINT  'ENTER  Y  AXIS  LABEL'NINPUT  LY$ 

430  !  CALCULATE  THE  RESOLUTION  FOR  ThE  PLOT 

440  !  RX  IS  THE  RESOLUTION  FOH  THE  X  AXIS 
450  !  rY  IS  THE  RESOLUTION  For  THE  Y  AXIS 
460  KS«I0000 

4/0  IP  AdS( XX)>ABS(XM)  THEN  RX-ABSC XX)  ELSE  RX-ABS(XM) 

4 BO  HX=KS/RX 

490  IP  ABS( YY )>ABS( YM )  THEN  RY-ABSCYY)  ELSE  KY=ABS( YM ) 

500  RY=RS/kY 

520  !  START  PLOTTING 

530  0X=I05*  !  OUTPUT  AND  INPUT  IEEE  ADDRESS 

540  I%=I05* 

550  IP  PC%=  U  GOTO  600 

555  PRINT  CHRS( 7*) I'PEN  CHANGE?  (Press  screen  to  continue )'\WAIT  FOh  KEY 

556  K*=KEY 
5/0  ! 

600  PRINT  alOX.'INISIl'  !  INITIALIZE  PLOTTER 

710  ! 

720  !  SET  LIMITS  OP  THE  PLOT 

730  PX  =  1 000* ( XL+2 )\PY= I 000* <  YL+ .5)  !  CALCULATE  PLOTTER  LIMITS 

740  PRINT  WOX.'IP  2000, 500, '«PX»'.'»PY;'«'  !  SEND  PLOTTER  LIMITS 

745  IF  PCX<>  IX  GOTO  900  !  SKIP  X  AXIS  FOR  MULTI-PLOT 

760  !  PUT  UN  X  AXIS 

770  PRINT  60*,'PA2000,500lPD»XTl' 

730  FOR  1=1  TO  XL 

790  PHINT  «0X , ' P  H I  000 ,  0 1  XT  I  ' 

300  NEXT  I 

dio  print  eox.'pui' 

320  PRINT  «K)X,'PA4000,  I OOlLB' »LXS»CHRSC3 ) »' »'  !  WRITE  LABEL  ON  X  AXIS 

830  DX=( XX-XMJ/XL 

840  print  uox.'pai 400,3001' 

850  FOR  1=0  TO  XL  !  SCALE  NUMBERS 

800  Z=*XM*DX*I 

801  ZS=NUM$(Z,'S**>#*.##') 

3/0  PRINT  «K)X, 'LU'lZSJCHRS(3)l'«'»'PKI  18,01' 

830  NEXT  I 

890  !  END  OF  X  AXIS  NOW  PUT  ON  Y  AXIS 
900  XVX=2500X-PC4*5004 

902  print  mox.'pa'.xvx.'.sooipoiyti' 

910  FOR  1=1  TO  YL 

920  PRINT  30X, 'PRO, 1 0001 YT I ' 

930  NEXT  I 

940  PRINT  hOX.'PUIUIO,  HPn-300,',-YL*750,'lLB'»IYXl'=  '  I LYS  ICHRS  ( 3  )  I  '  «  ' 

950  UY=( YY-YM1/YL 

900  PRINT  «OX,'PA', (XVX-I50X),',50»' 

970  PUR  1=0  TO  YL 

980  Z=Ym+DY*I 

901  ZS=NU«S(  Z,  'S»»#.  **#'  ) 

990  PRINT  ttOX,'Lu'»ZS»CrtH$(3>l'»PR0,  1G0»' 
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1000  NEXT  I 

1020  !  END  OP  Y  AXIS 

1022  IF  PCSSOI  GOTO  1090 

1030  !  WHITE  LABEL  ON  AXIS 

1040  FOW  I L=2  TO  4 

1050  PL=PY+I000-(IL-2>*250 

1060  PHINT  «O%,'PUIDII,0lPA2G00,'IPLl'lLB'lB2$CIL>lCHrtS(3)l'l' 
1070  NEXT  IL 
10/1  PL=PY*250 

1072  PHINT  «0%, 'PUIDI I ,0lPA2G00, ' I PLl ' ILB' 1PHS  ICHKS <  3) I' »' 

1090  XM3!-INT(XM*HX)SXX5;=INT(XX*mX) 

1092  YM4=INT(Ym*HY)\YY4=INT(YY*HY) 

1095  PHINT  «0%,'SC'tXM%l-','lXX%l‘','lYrt%«','IYY4»'l'!  SCALE  PLOT 
1100  !  PLOT  THE  DATA 

1124  1=0 

1125  PHINT  «04,'PA'IINT<X<I)*HX)l','IINTCY(I>*KY>l'lPD' 

1128  FOH  1  =  1*  TO  NP%-I% 

1130  PkINT  •0*,'PA'IINT<X<I>*HX)I','«INT(Y( I)*HY)I'IPD' 

1140!  IF  I  <>  1  THEM  PHINT  «0%,X  ( I  )*KX» ','  I  Y(  I  >*HY 
1150!  IF  Iol  AND  I<NP4  THEN  PHINT  4104, ',' 

1180  NEXT  I 

1190  PHINT  CHHS<27)|'[ IOD'ICHhS<27)I'IIK'I 

1200  PHINT  <*0%,-'  l PUI PH  I  00,0101  I ,  OtLb'  1 1  Y4ICHH  S  (3  ) I' I  INI ' 

1210  PC%=PC%+U\PHINT  'PLOT  DONE!  NEXT  CHANNEL' 

1220  PHINT  «0*,'PU»' 

1230  GOTO  100 
1240  PHINT  <10%,  '  INI' 

1260  CLOSE  I 

1270  IF  EHH  <>  0  THEM  PHINT  'EHH=' IEHH, ' EHL=' lEHL 
1280  END 
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